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Plasticity for the Aerodynamicist 


WILLIAM PRAGER* 


Brown University 


ABSTRACT 

With the increasing need for exact methods of stress analysis 
beyond the elastic range, the interest in the mathematical theory 
of plasticity is rising steadily. The present paper contains a 
survey of the classical theory of plasticity (Saint Venant-Mises 
theory) as applied to problems of plane plastic flow. The pres- 
entation is kept as elementary as possible, the only prerequisites 
being some knowledge of strength of materials and the rudi- 
ments of the theory of elasticity. The most important classical 
results, as well as some modern developments, are presented 
with particular emphasis on the analogies between plasticity 
and fluid dynamics. 


(1) InrRopucTION 


_ MATHEMATICAL THEORY OF PLASTICITY owes its 
development to the demand for more realistic 
methods of determining safety factors of structural 
members and machine parts and to the need for better 
control in technological forming processes such as 
rolling, drawing, or extruding. The plastic strains 
encountered in these two fields of application are of 
different orders of magnitude. In a hollow sphere 
under gradually increasing interior pressure, for in- 
stance, the plastic zone begins to form at the interior 
surface and gains in depth with further increase in 
pressure. Until the plastic zone penetrates to the ex- 
terior surface, however, it remains contained in an 
elastic shell, and the plastic strains remain of the 
same order of magnitude as the elastic strains. Large 
plastic strains become possible only after the plastic 
zone has broken through to the exterior surface. The 
use of the sphere as pressure container is obviously re- 
stricted to the phase of contained plastic deformation. 
On the other hand, the blowing of the sphere into a 
form is possible only in the phase of unrestricted plastic 
flow. In problems of unrestricted plastic flow the re- 
versible elastic deformations can usually be neglected 
in comparison with the permanent plastic deformations; 
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in problems of contained plastic deformation, however, 
this is not admissible, as a rule. 

Problems of contained plastic deformation seemed 
to call for an extension of the theory of elasticity; 
problems of unrestricted plastic flow, for an extension 
of fluid dynamics. Accordingly, different theories of 
plasticity were developed for these two fields of ap- 
plication. The theory of plastic deformation is based on 
a relation between stress and plastic strain; the theory 
of plastic flow, on a relation between stress and rate of 
plastic strain. Both theories assume that plastic de- 
formations can occur only after the yield limit has been 
reached—i.e., after the intensity of stress has reached a 
certain critical value. Numerous measures for the 
intensity of stress have been proposed; in the follow- 
ing, the maximum shearing stress will be used for this 
purpose (Tresca'). Once the yield limit has been 
reached, continued plastic flow may require continual 
increase of the intensity of stress (work-hardening 
materials) or may take place under constant intensity 
of stress (perfectly plastic materials). In this introduc- 
tion to the mathematical theory of plasticity, only 
perfectly plastic materials will be considered. More- 
over, the discussion will be restricted to the theory of 
plastic flow (theory of Saint Venant and Mises?) which 
is of greater methodological interest to the aerody- 
namicist. 


(2) Basic EQuaTIons 


The examples discussed in the following will all be 
concerned with quasi-static, two-dimensional flows of 
incompressible, perfectly plastic materials. Let oz, oy, 
o, denote the normal stresses and Ty, Tz, Ty: the tan- 
gential stresses with respect to the rectangular Car- 
tesian coordinates x, y, z. If the z axis is chosen normal 
to the plane of flow, 


= (oz + 
(see for instance, reference 3, p. 183). Thus, the state 


Tz: = Ty = 0 
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of stress at a point is completely defined by the values 
that o,, o,, and 7,, assume at this point. 

Inertia forces can be disregarded in a quasi-static 
flow. In the absence of other body forces, the stresses 
will then satisfy the equilibrium conditions 


Oren 


Moreover, the stresses will satisfy the yield condition 
(oz — oy)? + = 4k? (2) 


which expresses the fact that the maximum shearing 
stress equals the yield stress k in pure shear. 

Under appropriate boundary conditions, Eqs. (1) 
and (2) will be sufficient for the determination of the 
stresses. It is with these statically determinate problems 
of plasticity that the following discussion is concerned. 

For the following, it will be useful to discuss the 
stresses transmitted across surface elements that are 
perpendicular to the plane of flow. Fig. 1 indicates the 
notations and sign conventions adopted here, the 
stresses shown being those transmitted from the un- 
shaded to the shaded side of the surface element. The 
variation of NV, and T,, with a is represented by Mohr’s 
circle (Fig. 2a). The points X and Y correspond to a = 
0 and a = w/2 and have the coordinates o,, —7,, and 
Sy, Try, Tespectively. The segment XY is a diameter of 
the circle; according to the yield condition (2), the 
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radius equals k. The abscissa of the center represents 
the mean normal stress 


o = 1/2 (or + oy) (3) 


The points M and N represent the stresses transmitted 
across the principal planes of stress, and the points | 
and II represent the stresses transmitted across the 
planes of maximum shearing stress. The orientations 
of these planes and stresses are shown in Fig. 2b, 
The plane of maximum normal stress (a = a;) will be 
called the first principal plane; the plane of minimum 
normal stress [a = a + (2/2)], the second principal 
plane. When these planes are rotated in the counter- 
clockwise sense through 45° about their line of inter- 
section, they are made to coincide with the first and 
second shear plane, respectively [a = a; + (1/4) = 86 
and a = a, + (82/4) = 6+ (x/2)]. The traces of the 
first principal plane and the first shear plane in the 
x, y plane are called the first principal direction and the 
first shear direction, respectively. The first shear lines 
are defined as those curves in the x, y plane which have 
everywhere the first shear direction; the second shear 
lines are their orthogonal trajectories. 
With w = o/2k, it follows from Fig. 2a that 


o, = k(2w + sin 26) 
oy = k(2w — sin 26) 
Ty = —k cos 20 


(4) 


When these expressions, which satisfy the yield condi- 
tion, are substituted into the equations of equilib- 
rium, the following equations are obtained: 


w, + 6, cos 26 + 6, sin 20 = “i 


w, + 6, sin 20 — 6, cos 20 = 0 (5) 


Here, the subscripts x and y indicate differentiation 
with respect to these variables. 

To arrive at a physical interpretation of Eqs. (5), take 
the axes of x and y in the first and second shear direc- 
tions at a generic point P. The angle @ at this point 


Fic. 2a (left). 


Fic. 2b (right). 


then a 
respec 
first a1 
tion in 
respect 
the fol 


(¢ 


Thus, 
wt 


(3) 
Whe 


a shea 
followi 
lines it 
series | 
Since ¢ 
second 
is the 
shear | 
line at 
approx 
shear | 
B. If 
which 
w — 6 
wo + 6 


approx 


From t 


Simila1 
the val 
6 at 1, 
as the 

This 
variou: 
in the 
given 1 


This re 
first sk 
secant 
may b 


Anoth 
by an: 
of the 
gent tc 


| 
1 
Ox Oy oy Ox (1) 
: 
I 
k 
N 
N M No Y 
x 
X 
g=2kw 
°x 


PLASTICITY FOR THE AERODYNAMICIST 255 


then assumes the value 2/2, and differentiation with 
respect to x and y is equivalent to differentiation in the 
first and second shear directions at P. If differentia- 
tion in these directions is denoted by d/ds,; and d/ds2, 
respectively, Eqs. (5), written for the point P, assume 
the following form: 


(d/dsi)(# — 6) = 0,  (d/ds2)(w + 6) = 0 (6) 


Thus, # — @ is constant along the first shear lines, and 
+ is constant along the second shear lines. 


(3) APPROXIMATE CONSTRUCTION OF SHEAR LINES 


When w and @ are given along a curve c, which is not 
a shear line (Cauchy problem), Eqs. (6) suggest the 
following approximate determination of the shear 
lines in the neighborhood of c. Let A, B, C,...bea 
series of not too widely spaced points on ¢ (Fig. 3). 
Since @ is given along c, the tangents to the first and 
second shear lines at A, B, C,... may be drawn. If 1 
is the point of intersection of the tangent of the first 
shear line at A and the tangent of the second shear 
line at B, the point 1 may be considered as a first 
approximation to the point of intersection of the first 
shear line through A and the second shear line through 
B. If subscripts are used to indicate the points at 
which w and @ are to be evaluated, the constancy of 
w — 6 along first shear lines and the constancy of 
w + @ along second shear lines furnish the following 
approximate equations: 


+ = we + 


— 0, = — 
From these, w; and 6; can be determined: 


= + wp — O04 + (7) 
0; = 1/2 (04 + On — wa + wp) 


Similarly, we and 42 (see Fig. 3) can be determined from 
the values of w and @at Band C. The values of w and 
6 at 1, 2, ... can then be used in the same manner 
as the given values along ¢, etc. 

This somewhat crude procedure can be improved in 
various ways. For instance, instead of being drawn 
in the direction defined by @,, the segment Al may be 
given the direction defined by 


(04 + 61) = (30, + 05 — wa + wp) 


This refinement consists in replacing the tangent of the 
first shear line at A by a better approximation to the 
secant Al. In a similar manner, the segment Bl 
may be drawn normal to the direction defined by 


(Oe + 6:1) = '/4 (04 + 303 — wa + wp) 


Another refinement consists in replacing the segment A1 
by an arc of an approximation to the circle of curvature 
of the first shear line through A. This circle is tan- 
gent to the first shear line at A and has the radius 


Al _ 24B|cos (6 — a) | (8) 


Op — 94 + we — 


6; — 04 


Fic. 3. 


Fic. 4. 


(see Fig. 3). If the right-hand side of Eq. (8) is posi- 
tive the center of the circle lies to the left of an ob- 
server who proceeds along the first shear line through A 
in the direction of increasing values of w + 6. Similarly, 
the segment B1 may be replaced by an arc of the circle 
that is tangent to the second shear line at B and has the 
radius 


2AB\sin(@ — a)| (9) 


= % — 04 — wet ws 


If the right-hand side of Eq. (9) is positive, the center 
of this circle lies to the right of an observer who 
proceeds along the second shear line through B in the 
direction of increasing values of w — @. 

When the curve c itself is a shear line, these methods 
of constructing the net of shear lines in the neighbor- 
hood of ¢ fail. Indeed, if ¢ is a first shear line, for in- 
stance, the segment A 1 is a segment of c, and the point 1 
coincides with B. To define the net of shear lines in 
the neighborhood of a first shear line c, a second shear 
line d, which intersects c, must be given (Riemann prob- 
lem). Along the first shear line c, the quantity 6 — w 
has a constant value. Since the variation of @ along c 
is known from the shape of this curve, the variation 
of w along c is known, too, if the value of w at one point 
of ¢ is known. A similar statement applies to the 
second shear line d. If 0 is the point of intersection of 
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c and d (Fig. 4) and if A and B are neighboring points 
on d and ¢, respectively, the point 1 can be constructed 
from the data at A and B, and the values of , and 4, 
can be determined in precisely the same manner as 
above. Application of the same procedure to 1 and 2 
(Fig. 4) leads to point I, etc. 

It follows from the preceding discussion that the 
shear lines are the characteristics of the differential Eqs. 
(5). They play the same role in these problems of 
two-dimensional plastic flow as the Mach lines in 
problems of two-dimensional supersonic flow: _Dis- 
continuities in the derivatives of w or @ and, hence, 
discontinuities in the derivatives of the stress compo- 
nents can occur only along shear lines; if the stress dis- 
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tribution in two adjacent regions is described by dif. 
ferent analytic expressions, the boundary between these 
regions is a shear line. 


(4) BounDARY CONDITIONS 


In Fig. 5, let the curve Os represent part of the 
contour of a two-dimensional domain, which lies on 
the shaded side of this curve. Take O as the origin 
of the arc length s and choose the sense of increasing s 
so that the domain remains to the left of an observer 
who moves in this sense along the contour. The nor- 
mal and shearing stresses transmitted across the con- 
tour will be denoted by N and T, respectively. The 
stress NV will be counted as positive if it constitutes a 
tension exerted on the domain under consideration, 
and the stress 7, if it has the direction of increasing s, 
From a comparison of the sign conventions indicated in 
Figs. 1 and 5 and from the Mohr diagram of Fig. 6, it is 
seen that 


N/k = w+ sin 2(6 — a) (10) 
T/k = cos 2(6 — a) (11) 


Thus, if NV and 7 are given at a point of the boundary 

where the exterior normal makes the angle a with the 

positive x axis, the values of 6 and w at this point are 
a + '/, are cos (T/k) 


given by 


Because of the ambiguity of sign in the first of these 
equations, there are two pairs of values 6, w for each 
given pair of values 7, N. This ambiguity springs 
from the quadratic character of the yield condition. 
Indeed, if N’ denotes the normal stress transmitted 
across a surface element that is normal to the contour 
at P, the yield condition (2) is equivalent to 


(N — N’)? + 47? = 4k? 


When WN and T are given, this equation furnishes only 
the absolute value of the difference N — N’ but not its 
sign. Which of the two possible values of NV’ must be 
taken in a specific case cannot be decided from the 
local conditions at a boundary point but depends 
on the complete boundary conditions. Fig. 7 shows 
a particularly simple example. It is not possible to 
decide from the local conditions at P whether the 
normal stress N’ transmitted across a surface element 
that is normal to the contour at P equals p + 2k or p — 
2k. From the complete boundary conditions of the 
problem it follows, however, that V’ = p — 2k. 


ll 


(13) 


(5) TyprcaL PATTERNS OF SHEAR LINES 


The following remarks concerning straight shear lines 
often prove useful in the construction of shear line 
patterns. Along a straight shear line, w, as well as 8, is 
constant, and any straight shear line is embedded in a 
family of straight shear lines. Indeed, along a first 
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shear line, for instance, @ — wisconstant. If this shear 
line is straight, @ is constant along it. Thus, w is also 
constant along a straight shear line. Moreover, 
since the state of stress is completely defined by the 
values of w and 6, the same state of stress prevails in 
all points of a straight shear line. 

To prove the second statement concerning straight 
shear lines, consider a straight first shear line AB 
(Fig. 8), a neighboring first shear line CD, and the 
second shear lines AC and BD. Since AC and BD 
are second shear lines, 


wa + O4 = we + %, 


Since AB is a straight shear line, 


+ 03 =wp+ (14) 


wa = Wz, 04 = 62 (15) 
From Eqs. (14) and (15), it follows that 
we + O¢ = wp + O (16) 
Since CD is a first shear line, 
we — = wp — (17) 
Thus, 
= Wp, 9c = (18) 


Since these relations are valid for any two points, C 
and D, on this neighboring shear line, this line, too, 
must be straight. Moreover, Eqs. (14), (15), and (18) 
show that w + @ is constant throughout this field of 
straight first shear lines. Similarly, w — @ is constant 
in a field of straight second shear lines. 

In a field of parallel straight first shear lines, @ is con- 
stant in addition to w + 6. Accordingly, w, as well as 
#, is constant throughout the field, which, therefore, 
is a field of constant stress. 

If the state of stress is constant along an arc of a 
curve, this arc lies in a field of constant stress. (A 
straight segment along which the state of stress is con- 
stant might be a segment of a straight shear line and, 
hence, not lie in a field of constant stress.) 

In the following, a region where only one family of 
shear lines consists of straight lines will be called a fan; 
if, in particular, these straight lines intersect in the 
same point, the term centered fan will be used. Fans 
obviously play the same role in plasticity as simple 
waves in supersonic flow. 

Fig. 9 illustrates a problem in plane plastic strain 
which can be solved in terms of fans and regions of 
constant state. The blunt wedge FABG is subjected 
to uniformly distributed pressure along AB. Along 
each of the straight segments FA, AB, and BG, the 
state of stress is constant. Moreover, these segments 
are free from shearing stresses; they therefore con- 
stitute lines of principal stress which the shear lines 
intersect under +45°. Accordingly, each of these seg- 
ments lies in a region of constant state which has the 
shape of an isosceles right triangle (Fig. 9). The 


transition between the regions of constant state ADF 


Fic. 8. 


and ABC is provided by a fan centered at A, and the 
transition between the regions of constant state BCA 
and BGE is provided by a fan centered at B. The 
material below the line FDCEG may be treated as a 
rigid body. 

To decide which of the shear lines thus determined are 
first and second shear lines, consider the deformation 
of the shaded square mesh in Fig. 9. The lower side of 
this mesh forms part of the boundary of the rigid re- 
gion and therefore does not move at all, while the upper 
side moves to the left. The diagonal that is parallel to 
AF is shortened and the other diagonal is lengthened 
during this deformation. Since the principal stress 
perpendicular to AF vanishes, the principal stress 
parallel to AF must thus be negative (compression). 
From this it follows that the first shear lines in the 
region ADF are parallel to DF. The directions of the 
first shear lines in the other zones are then determined 
by continuity considerations. In Fig. 9, the first and 
second shear lines are indicated by full and dotted lines, 
respectively. 

The intensity of pressure necessary to produce yield 
ing in accordance with the shear line pattern of Fig. 9 
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can be found as follows. With the notations intro- 
duced in Section 4, the boundary conditions along AF 
are: 


N=T=0 
The compressive stress N’ which is required to fulfill 
the yield condition is therefore given by N’ = —2k. 


Thus, 
o= (V+ N’)/2 
2k 


along AF and, hence, in the region ADF. In the fan of 
second shear lines, ACD, the quantity w — @ is constant, 
while @ increases by f in the transition from AD to AC. 
Accordingly, w also must increase by £ in this transition. 
Thus, 


2) 


w=B-—'/, O= 4/4 


in the region ABC and, hence, along AB. With these 
values, the second Eq. (4) gives 


o, = k(2w — sin 26) = —2k(1 + 8) 


This is the pressure necessary to produce yielding. 
Fig. 10 shows the pattern of shear lines in a wedge 
that is drawn through a lubricated die. In the re- 
gion ABGECDF the pattern is the same as in Fig. 9. 
From the values of w and @ along the shear lines CD and 
CE the pattern in the region CDHE can be determined 
(Riemann problem). The fans DFIH and EGJH 
are formed by the normals of the arcs DH and EH; 
they are cut off by the orthogonal trajectories FI and 
GJ. The sequence HIM/J is a region of constant state. 
The shear lines in the regions FIK and G/L are deter- 
mined from the values of w and @ along the shear lines 
FI and GJ and the constant values of @ along FK and 


(19) 


GL (mixed problem). The fans JKNM and JLOM 
are formed by the normals of the arcs JK and JL; 
they are cut off by the orthogonal trajectories MN 
and MO. The triangles KNP and LOQ, finally, are 
regions of constant state. Sokolovsky (reference 4, 
p. 187) gives the pattern of shear lines for a nonlubri- 
cated die, assuming that the intensity of the frictional 
shearing stress has a given constant value. 


(6) Lriuitinc Lines AND LINES oF DISCONTINUITY 


When the value of w and @ are given along an arc AB, 
which is not a shear line, the methods developed in 
Section 3 will furnish the shear lines within a certain 
region that includes this arc. As a rule, this region is 
the curvilinear quadrangle ACBD bounded by the 
shear lines through A and B (Fig. lla). If, say, the 
first shear lines have an envelope, it is possible, how- 
ever, that the first shear line through A meets this 
envelope before intersecting the second shear line 
through B. In Fig. 11b, for instance, the arc AB is 
assumed to lie in a noncentered fan of first shear lines. 
The envelope of these straight shear lines is denoted by 
e, and the first shear line through A is tangent to this 
envelope at E, while the second shear line through B 
meets this envelope at D. In this case the data along 
AB determine w and @ and, hence, the stresses in the 
region AEDBC. 

With reference to Figs. 1la and 11b, now let the values 
of w and 6 be given along an arc A’B’ that contains the 
arc AB. In the case of Fig. lla, the stresses are then 
determined in the region A’C’B’D’, which contains the 
region ACBD in its interior. In the case of Fig. 11b, 
however, the stresses are determined in the region 
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Fic. lla (left). 


A'E’D’B'C’, which has the boundary are ED in 
common with the region AEDBC. The envelope of a 
family of shear lines thus plays the role of a limiting 
line beyond which the stress analysis cannot be con- 
tinued. 

In the case illustrated by Fig. 11b, any second shear 
line (orthogonal trajectory of the straight first shear 
lines) has a cusp at the point where it meets the en- 
velope of the first shear lines. (It can be shown that, 
generally, the envelope of one family of shear lines 
is the locus of the cusps of the shear lines of the other 
family.) Thus, the shear lines of neither family can 
be extended beyond the envelope of the shear lines of 
one family. 

In Fig. 12, let AB and CD be two neighboring lines 
of a noncentered fan of first shear lines with the en- 
velopee. Now, the quantity w + 6 is constant through- 
out such a fan, and w, as well as 6, is constant along any 
straight shear line. Since @ increases by Ag (Fig. 12) 
in the transition from AB to CD, the quantity w must 
decrease by Ay. At the point A, the derivative of w 
or 6 in the direction normal to AB has, therefore, the 
absolute value (AB)-!. This shows that in the points 
of the envelope e, the derivatives of w and @ and, 
hence, the derivatives of the stress components in the 
direction normal to the envelope become infinite. (It 
can be shown that this statement also holds for the 
envelope of a family of curved shear lines.) 

The equations of equilibrium which involve the first 
derivatives of the stress components become meaning- 
less in the points of a limiting line. The shear line 
pattern constructed from the boundary data must 
therefore be broken off and replaced by another solution 
before it reaches a limiting line. This, obviously, can- 
not be done in a continuous fashion but requires dis- 
continuities in some stress components. In Fig. 13, 
let a-a be a line of discontinuity. Equilibrium re- 
guires the stress components N and T to be continuous 


Fic. 11b (right). 


across this line, but the stress component N’ need not 
be continuous. Since the yield condition (13) is to be 
fulfilled on both sides of the line of discontinuity it is 
easily seen that 


— = +4V — T? (20) 


Fic, 12. 
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A well-known example of such a line of discontinuity 
is afforded by pure flexure in the plastic range. The 
distribution of bending stresses over the height of the 
beam is indicated in Fig. 14. The line a-a is a line of 
discontinuity; the stress NV’ equals the yield stress in 
tension on one side of this line and the yield stress in 
compression on the other; the stresses N and T vanish. 

A less trivial example has recently been discussed by 
the present author. Fig. 15a shows the shear line 
pattern in an infinite wedge that is exposed to uniformly 
distributed pressure on one of its sides. For obtuse 
wedge angles the pattern is the same as the pattern 
in the neighborhood of the point B in Fig. 9. For acute 
wedge angles, however, the regions of constant state 
represented by the isosceles right triangles ABD and 
BEC in Fig. 15a would overlap and the shear line ADEC 
would double back on itself (Fig. 15b). This can be 
avoided by the discontinuous solution shown in Fig. 
15c. The regions of constant state adjacent to AB 
and BC meet in a discontinuous manner along the bi- 
sector of the wedge angle BD. The pressure that 
produces yielding in accordance with this pattern is 
found to be 

oy = —2k(1 — sin (21) 


(see reference 4). For obtuse wedge angles the pres- 
sure is given by Eq. (19). For a wedge angle of 90° 
(8 = 0) both formulas give ¢, = —2k. 


(7) VeELocity Frecps 


The theory of plane plastic flow assumes that the 
shear lines indicate not only the directions of maximum 
shearing stress but also the directions of maximum rate 
of shear. This means that, at a generic point, the 
unit extensions in the first and second shear directions 
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must have equal values. The condition of incompres- 
sibility, on the other hand, requires that the sum of 
these unit extensions vanishes. Accordingly, each of 
these two unit extensions must vanish separately, 
The relations expressing this fact are the differentia] 
equations for the determination of the velocity field, 
These equations are most readily established in the 
case of a field of constant stress. In Fig. 16 the co- 
ordinate axes are taken parallel to the straight shear 
lines of such a field, and the velocity components in the 
first and second shear directions are denoted by U and 
V, respectively. The conditions that the unit exten- 
sions in the shear directions vanish are 


=0, dIV/dw=0 (22) 


Thus U is constant along any first shear line and V along 
any second shear line. 

If, in a field of constant stress, the velocity vector is 
given along an arc AB that intersects no shear line 
more than once, the velocity field is obviously deter- 
mined throughout the rectangle ACBD formed by the 
shear lines through A and B (Fig. 16). The velocity 
vector at C, for instance, has the same U component 
as the velocity vector at A and the same V component 
as the vector velocity vector at B. (If the arc AB 
intersects the same shear lines repeatedly, the velocity 
vectors cannot be given arbitrarily along AB.) This 
construction of the velocity field shows that the shear 
lines are the characteristics of the differential equations 
for the velocities, as well as of those for the stresses. 

To establish the differential equations for the velocity 
field in the general case of curved shear lines, consider 
two neighboring points, A and B, of a first shear line 
and the second shear lines through A and B. Further- 
more, defining the positive directions along these shear 
lines as the directions of increasing values of w + 0 
and w — 6, respectively, denote the velocity components 
in these directions by U and V. If V = 0, the unit 
extension in the direction of the first shear line is ob- 
tained as the derivative dU/ds, of U with respect to 
the arc length of the first shear line (Fig. 17). On the 
other hand, if U = 0, the unit extension in the direc- 
tion of the first shear line is found to be — V/R:, where 
R, is the radius of curvature of the first shear line, con- 


Cc 


Fic. 15a (left). Fic. 15b (center). Fic. 15c (right). 
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sidered positive when the center of curvature lies to 
the left of an observer who proceeds along this line in 
the direction of increasing values of w + 6. Similar 
considerations apply to the unit extension in the direc- 
tion of the second shear line. The conditions for the 
vanishing of the unit extensions in the directions of the 
shear lines therefore are 


(dU/ds:) — (V/R:) = 0, (dV/dsz) — (U/R:) = 0 
(23) 


These equations show that U is constant along any 
straight first shear line (Ri = ©) and that V is constant 
along any straight second shear line. Also, if V van- 
ishes along a first shear line, then U is constant along 
this line. Similarly, if U vanishes along a second shear 
line, V is constant along this line. 

If U and V are given along an arc AB that is not a 
shear line and does not intersect any shear line more 
than once, the velocity field inside the curvilinear rec- 
tangle bounded by the shear lines through A and B can 
be determined by means of Eqs. (23). In particular, if 
U = V = 0 along the arc AB, then U = V=0 
throughout this rectangle. The boundary of a region 
that remains at rest must therefore be formed by shear 
lines. 


(8) Prastic BouNDARY LAYERS 


The theory reviewed in the preceding sections re- 
sembles the theory of perfect fluids insofar as it neglects 
viscosity effects. The equations of motion of a perfect 
fluid are obtained from the equations of motion of a 
viscous fluid by simply setting the viscosity » equal to 
0. According to Prandtl,* physically important effects 
are lost by this outright procedure, which should be re- 
placed by the more subtle transition to the limit, 


As the present author pointed out (reference 7, pp. 
12 and 28), similar considerations may be important 
in plasticity. According to this point of view, the 
Saint Venant-Mises theory discussed in the preceding 
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should be considered as a limiting case of the Bingham 
theory. According to the first of these theories, a 
specimen in pure shear yields under a constant shearing 
stress k, the intensity of which is independent of the 
rate of shear. According to the second theory, the 
yield stress in pure shear exceeds k by an amount pro- 
portional to the rate of shear.* For more general 
states of stress, the relations between stresses and strain 
rates were established by Hohenemser and Prager.® 
For an incompressible material in plane flow, the 
mean pressure p and the maximum rate of shear IT are 
defined by 


—"/2 + oy) 


where u and v are the velocity components in the direc- 
tions of the coordinate axes. The relations of Hohenem- 
ser and Prager then yield the following equations for the 
plane flow of a Bingham material: 


“SyL ox 
(24) 


where uz is the coefficient of viscosity. For k = 0, these 
equations reduce to the well-known equations for the 
slow plane flow of a viscous fluid (reference 10, p. 22); 
for » = 0, on the other hand, they reduce to a less 
well-known form of the equations for the plane flow of a 
Saint Venant-Mises material. 

The transition to the limit, u — 0, leads to the equa- 
tions of the Saint Venant-Mises theory in the “free 
steam” and to “boundary-layer equations’ in the 
neighborhood of solid boundaries. For the case of a 
straight boundary along the x axis, the following bound- 
ary-layer equations were recently established by Old- 
royd:!! 


b + 2k[(Ou/dx)/(Ou/dy)] = f(x) 


(25) 
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Oldroyd has succeeded in integrating these equations in 
several cases where the material outside the boundary 
layer remains rigid. No solutions of Eqs. (24) seem 
to be known as yet in the more interesting cases where 
the boundary layer borders on a region of plane plastic 
flow of the Saint Venant- Mises type. 
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Letters to 


Dear Sir: 

This refers to the paper by E. H. Schuette, entitled “Buckling 
of Curved Sheet in Compression and Its Relation to the Secant 
Modulus,’”’ published in the January, 1948, JourNaL. Mr. 
Schuette’s results are interesting in their indication that buckling 
occurs at a relatively fixed value of critical strain, whether in the 
elastic range or the inelastic range. 

The theoretical buckling load of a cylindrical shell axially 
loaded, under the usual assumptions, gives for the strain at which 
buckling occurs in the elastic range, the formula 


écr = 0.60t/R 


for the value of Poisson’s ratio of about 0.3. 

The results indicated by Mr. Schuette range from about 51 
to 135 per cent of this theoretical value with the major part of the 
range, as shown by the dashed-line curve in Fig. 2, from 56 to 
100 per cent. The average value quoted in his conclusions and 
shown by the full line on Fig. 2 is about 70 per cent of the so-called 
theoretical value. 

In 1933, the writer, with W. M. Wilson, reported results of tests 
on thin cylindrical steel shells loaded axially as columns (‘‘The 
Strength of Thin Cylindrical Shells as Columns,” Bulletin No. 
255, University of Illinois Engineering Experiment Station, 
1933). The report summarized the results of tests on small 
carefully machined specimens, on large carefully machined speci- 
mens up to 14 in. in diameter, and on extremely large fabricated 
specimens with riveted joints and with welded joints. The 
riveted specimens ranged in size up to 80 in. in diameter by 6 ft. 
in height, and the welded specimens were of similar dimensions, 
but one series of welded specimens were 34 in. in diameter by 35 
ft. in length. 

In the report, the writer suggested the use of the critical strain 
interpretation of buckling for both the elastic and plastic ranges. 
This, I believe, is the first suggestion in the literature of such an 
interpretation and is, of course, the same as a secant modulus. 

The writer’s test results led to values of critical strain ranging 
from 18 to 26 per cent of the theoretical value for fabricated 
specimeus and from 31 to 45 per cent of the theoretical value for 

the relatively accurately machined specimens. The highest 


the Editor 


ratio for any specimen in the writer’s tests was 57 per cent. 
These results are in general qualitative agreement with those 
reported by Mr. Schuette, but they indicate that a more conser- 
vative attitude should be adopted in these problems unless some 
specific information is available regarding the material and the 
degree of curvature of the shell. 


N. M. NEWMARK 
Research Professor of Structural Engineering 
College of Engineering 
University of Illinois 


Dear Sir: 

In the January, 1948, issue of the JouRNAL, Lloyd H. Sloan, 
in a Letter to the Editor, p. 40, commented on the writer’s paper 
published in the June, 1947, issue. 

The point in question is whether or not, as Mr. Sloan stated 
in his original article, an aircraft flying at constant air speed will 
undershoot from the original glide in a gradient wind. The writer 
has shown in his paper that this is not the case and that the glide 
is, in fact, extended. 

Mr. Sloan has quoted from AP 129, the R.A.F. flying manual. 
In the 1941 edition of this publication, Fig. 78 and paragraph 267 
both clearly state that the gain in rate of descent being discussed 
is associated with a loss of air speed, resulting from a sudden de- 
crease in head wind. This must be considered as a gust condition 
rather than one of steady wind gradient. 

With regard to pilots’ reports concerning this phenomenon, 
these must be treated with considerable caution. Without any in- 
strumental record or quantitative evidence, it is not possible to 
say through what form of wind structure the descent took place. 
In a gust condition, with loss of air speed, there is, pf course, a 

steepening of the glide. One R.C.A.F. pilot has reported to this 
writer that he observed an extension of his glide under what he 
took to be a condition of steady wind gradient. 


BERNARD ETKIN 
University of Toronto 


This 
of the 
promp' 
cipal s 
indust1 
advanc 
indeed 
at this 

Fact 
dents 
gear fe 
mercia 
combi! 
nantly 
in desi 
been i 
to sho 
are sel 
diagra 
aircral 
behav: 

In 
in apf 
heed 
explar 
strut 
deflec 
to the 
menta 
gears 
drag | 
occur! 

Also, 

skids 

tions | 
motio 

ing d 


q 
‘ 
I 
/ 
man 
1 
7/2 
henc 
weig 
are ¢ 
such 
cavi 
it. 
Pr 
Meet 
H. 
semt 


A Critical Study of Aircraft Landing Gears’ 


J. F. McBREARTYt 
Lockheed Aircraft Corporation 


ABSTRACT 


This report contains some of the results of a continuing study 
of the behavior of aircraft landing gears. This study was 
prompted by the general awareness that landing gears are a prin- 
cipal source of accidents and troubles throughout the aviation 
industry and by the desire to determine, if possible, technological 
advancement for future aircraft. Some of the findings to date, 
indeed, appear of sufficient importance to warrant presentation 
at this time. 

Factual data are presented showing that, exclusive of acci- 
dents due to crew error, weather, or unknown causes, landing- 
gear failures have been consistently responsible for more com- 
mercial and military accidents than all the rest of the structure 
combined, that landing-gear structural failures are predomi- 
nantly from drag loads, that this is in spite of steady increases 
in design loads, and that a large number of these failures have 
been in “normal” landings. Some statistical data are available 
to show that design conditions are amply severe and, in fact, 
are seldom realized perhaps even in “hard” landings. Strength 
diagrams show that landing-gear strengths on comparable 
aircraft are about equivalent and are all very high. Likewise, 
behavior characteristics are shown to be similar. 

In a critical search for an underlying cause of excessive loads 
in apparently normal landings with aircraft of all types, Lock- 
heed Constellation experience is utilized. Various possible 
explanations, such as the general strength level, fatigue, shock 
strut characteristics, extreme ground friction, and structural 
deflections, are considered and found to aggravate and contribute 
to the problem but fail to explain service failures. The funda- 
mental cause is deduced to be dynamic overloading of landing 
gears resulting from critical rates of application of wheel spin-up 
drag forces, which critical rates have a random probability of 
occurrence on any aircraft unless the motion is suitably damped. 
Also, at least one source of self-excited vibration arising from 
skids on wet runways is found to exist. Some possible solu- 
tions are discussed, and hydraulic damping of horizontal wheel 
motion is advocated. Typical test results show that such damp- 
ing does indeed offer a most effective solution. 


INTRODUCTION 


|= AIRPLANE LANDING GEAR is notoriously heavy, 
expensive, and troublesome and has been so for 
many years. Landing-gear weights generally average 
7'/. to 10 per cent of the airplane empty weight and, 
hence, offer a fruitful field for weight economy. This 
weight represents only the landing gear itself, and there 
are other indirect weight costs chargeable to the gear, 
such as the structural penalty required to provide a 
cavity to contain the gear, to cover it up, and to actuate 
it. Similarly, the dollar cost of landing gears has al- 


Presented at the Aircraft Design Session, Sixteenth Annual 
Meeting, I.A.S., New York, January 26-29, 1948. 

* The author wishes to acknowledge the assistance of Mr. K. 
H. Achterkirchen of the Lockheed Aircraft Corporation in as- 
sembling some of the statistics in this paper. 

Division Engineer. 


ways been large, involving much intricate machining 
and precision assembly. 

These prices would perhaps be tolerable if the 
landing gear gave a measure of service reasonably 
comparable to the dollars and pounds invested in it. 
But the reverse is true. There is a general and well- 
founded feeling throughout the aircraft industry that 
the landing gear is a continual source of trouble, exces- 
sive maintenance, and downright failure. 

Aware of the great potential that the landing gear 
offers in weight and manufacturing and operating cost, 
Lockheed has been conducting an informal study of 
landing-gear behavior for several years. It was soon 
realized that the problems involved required the 
efforts of engineering specialists on landing gears to 
correlate the experience available, and engineers were 
therefore assigned to such staff functions. Similarly, 
a decided lack of research on landing-gear and shock- 
strut behavior was soon apparent. Drop testing 
facilities were known to be available at the shock-strut 
manufacturers, but efforts to obtain research-type 
data from these sources failed because of pressure of 
their routine work and lack of funds. The Army has 
drop test facilities at Wright Field, and what research 
data are available have come from this source. How- 
ever, a considerable fraction of this work is routine 
product testing, and much too little research is possible 
on the Army equipment. The N.A.C.A. has no drop 
testing facilities, although it does have an elaborate 
impact basin for research work on water landings of 
hulls and floats. 

Therefore, Lockheed undertook to construct a drop 
testing machine in order to conduct experimental work 
on Lockheed landing gears and to permit the evaluation 
of some of the phenomena hitherto determined only in 
a qualitative manner. This machine (Fig. 1) is now 
completed and has many outstanding features, such 
as complete electronic instrumentation, simple simula- 
tion of wing aerodynamic lift, and great capacity 
With the help of additional funds from the Navy 
Bureau of Aeronautics, research already in progress 
will, it is hoped, make substantial contributions to the 
knowledge of landing-gear behavior. 

Although a great deal more work is contemplated 
for the future, this study has reached a stage such that 
a report covering some of the aspects is appropriate 
at this time, since the conclusions are believed to have 
considerable significance in the matter of structural 
failures and service troubles that continue to occur on 
all types of aircraft. 
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Fic. 1. Lockheed landing-gear drop testing machine. 


SERVICE STATISTICS 


While it is generally felt among aircraft engineers and 
operators that landing gears have a rather poor record 
in service, it is doubtful if many appreciate how poor 
this record is. Accident records and service reports 
are spread out over long periods of time, and usually 
the overall importance ‘of the collective data is not 
exploited. Contributing to this lack of attention to 
landing-gear service records is the fact that landing- 
gear failures rarely cause injury or loss of life and conse- 
quently most effort, in the U.S. at least, has been 
directed on other problems that more directly affect 
safety. However, landing-gear failures and troubles 
do result in enormous losses from direct damages, 
reduced revenue, and inoperative military aircraft, 
and therefore a critical review of these data is pre- 
sented, 

Statistical accident data available have been ex- 
amined, and the results are shown ia Figs. 2, 3, and 4. 
Civil aircraft accidents are classified by the C.A.B. 
into four cause factors: personnel errors, material 
(aircraft) failures, miscellaneous (weather, airport 
etc.), and undetermined. Material failures are further 
subclassified to four categories: power-plant failures, 
structural failures, handling qualities, and instruments. 
Handling qualities include all accidents attributable 
to inherent control characteristics and, together with 
instruments, constitute a small percentage of material 
failures. 


PERCENT OF MATERIAL FAILURES 


PERCENT OF MATERIAL FAILURES 
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DOMESTIC AIR CARRIER ACCIDENTS 
NOTES: 


1. MAIN GEAR INCLUDES RETRACTING MECHANISM. 
2. MATERIAL (AIRCRAFT) FAILURES INCLUDE POWER 
PLANT FAILURES, STRUCTURAL FAILURES, 
HANDLING QUALITIES, AND INSTRUMENTS. 
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AIRCRAFT LANDING GEARS 


TABLE 1 
Air Carrier. Some Landing Gear Accident Reports 
(Reference: CAA Air Carrier Operations by Extent of Injury and Damage, 1939 through 1945) 


Date Place Cause (Quoted) 

12/16/39 Caripito, Venezuela After a normal landing and while rolling on ground, left landing gear collapsed. 

8/15/40 Traverse City, Mich. After landing, landing-gear strut failed, ground-looped. 

6/5/41 Tampa, Fla. When in landing roll, right landing gear failed, caused by striking top of dyke on a 

i previous landing. 

11/11/41 Kansas City, Mo. When landing, right landing gear failed because of weakened condition of this member, 
causing a ground loop. 

9/17/42 Tallahassee, Fla. During the landing roll, the left landing-gear drag strut fitting failed. 

11/13/43 Kansas City, Mo. Pilot failed to level wings in cross-wind landing, and plane apparently was drifting at 
time of contact with runway. During the landing roll, the left landing gear collapsed. 

8/6/44 Watertown, S.D. Lower right drag strut failed through a faulty weld during a normal landing roll. 

7/15/45 Boston, Mass. In landing, the left landing gear drag strut failed, due probably to previous weakening. 


Structural failures, in turn, include the air frame and 
the complete landing gear with wheels, tires, brakes, 
etc. Fig. 2 isa plot of all material failures on domestic 
air-carrier airplanes from 1929 to 1944, the most recent 
year for which data are available. Complete landing- 
gear and air-frame failures are plotted in per cent of 
all material failures, but, for clarity, power plant, 
instruments, and handling qualities are not plotted, 
since they constitute the remainder. The curves 
show that the landing gear is second only to the power 
plant as an accident cause in the material failure 
category. Furthermore, while wheels, tires, and brakes 
are substantial contributors, the main gear without 
these items is responsible for more accidents than is the 
remainder of the entire airplane. In fact, in 1944 the 
main gear alone was responsible for a/J material failure 
accidents. The trend of these data also shows that 
the landing gear in general and the main gear in par- 
ticular have been responsible for an increasing fraction 
of accidents attributable to the airplane itself. 

The domestic nonair-carrier accident data shown in 
Fig. 3 reveal a similar situation, although power plant, 
handling qualities, and instruments constitute a larger 
share. The landing gear is responsible for twice the 
accidents due to the entire air frame, and the main gear 
alone far exceeds the air frame in every year of this 
15-year period. Incidentally, the retracting system 
cannot be charged with a large share of this trouble, 
since the data exclude retracting system failures from 
the main gear statistics. 

The available data on Army Air Forces’ accidents 
show similar features, as will be seen from Fig. 4. 
For the period 1943 to 1946, the landing gear was 
responsible for over three times the accidents due to 
the entire air fsame or about one-third of all failures 
ascribed to the aircraft itself—again second only to 
the power plant. 

The magnitude of the landing gear as a source of 
accidents is dramatically summarized in the following: 

“In this three months’ period, there were sixty 
cases of collapse of the landing gear. Thirty of 


these occurred during the normal landing roll.”’* 


* From: Supplemental Analysis No. 72, dated 3 July, 1945, 
A.A.F. Office of Flying Safety, Landing Gear Failures in Tactical 
Aircraft, Continental U.S., January through March, 1945. 


The underscoring in this quotation has been added by 
the author, and further discussion of this point is pre- 
sented later in this report 

The foregoing data cover only those cases that result 
in accidents and do not include the failures that have 
less serious consequences. Unsatisfactory reports show 
a similar indication, and it has been reported that, 
neglecting combat damage, trouble with the landing 
gear was responsible for more aircraft in unflyable 
condition during the war than all other causes com- 
bined except for the power plant, to which the landing 
gear was a close second. Air-line maintenance per- 
sonnel report similar experience, and landing-gear 
overhaul periods reflect this poor service record. 


SPECIFIC SERVICE REPORTS 


It is frequently claimed that the landing gear should 
be responsible for the bulk of structural failures and 
troubles because it takes more punishment, and it is in 
landing that the maximum skill of the pilot is required. 
It may be true that more skill is required to land an 
airplane than to fly it, but it remains to be proved that 
the pilot actually does a poorer job of landing than 
flying. Besides, obviously bad landings that produce 
structural failure are invariably classified as pilot error 
and are thus not included in the statistics presented 
on material failures. 

There might be grounds for the argument, however, 
if it were not for the fact that a large number of such 
failures occur in ‘“‘normal’’ landings. As will be pointed 
out later in this paper, the design landing conditions 
are so severe that, when a landing is made where the 
conditions even approach those used in design, the land- 
ing is described by any and all observers as bordering 
onacrash. Thus, for failures to occur in normal land- 
ings requires more than a casual explanation. They 
cannot be dismissed by minor discrepancies in work- 
manship, maintenance, or landing circumstances. 

Table 1 contains some of the examples found in the 
data available. These accidents to air-carrier airplanes 
reveal that it is not unusual for landing-gear failures 
to occur in “normal’’ landings. Table 2 is a similar 
list of examples found in Army Air Forces’ records. 

A study of service and accident records, exemplified 
in Tables 1 and 2, shows that failures seldom occur 
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TABLE 2 


Some A.A.F. Landing Gear Accident and Trouble Reports 
(Quoted from U.R. Digest) 


Date of 


U.R. Digest Model 
9/5/44 
Se 


10/20/44 B- 
12/5/44 AT- 
12/20/44 A- 
5/20/45 C- 
7/5/45 
7/20/45 B- 
8/5/45 C- 
8/20/46 C- 
12/5/46 P- 


Difficulty 
The left landing gear collapsed 
after a normal landing had 
been made. 


Right main landing gear, P/N 
65680, broke off as airplane 
touched the ground after a 
normal approach and landing. 
Investigation revealed the 
inner cylinder, P/N 65683, 
had sheared. 


Upon landing and rolling for 
several yards, the landing gear 
collapsed because of failure of 
trunnion bolt. 

During normal landing, right 
landing gear collapsed to the 
front. 


During landing, airplane ground 
looped as a result of failure of 
fitting assembly—landing gear 
strut attaching, P/N 5117828. 


A main landing gear strut, 


Bendix P/N 67270 (or 67271), 
collapsed. It was found that 
the pivot shaft and ring as- 
sembly, Bendix P/N 66584, 
was broken. 


Main landing gear collapsed 


upon landing. Gear broke 
just above flange where 
bayonet is bolted to strut. 


Bolt, landing gear upper attach- 


ing, P/N 1113280, failed dur- 
ing norma] landing. 


Airplane made a normal landing, 


and after approximately 2,500 
ft. of the landing roll had been 
completed, the left landing 
gear collapsed. 


While performing 50-hour in- 


spection, support-landing gear 
trunnion, P/N 117-33104, was 
found to be cracked through 
the trunnion upper bolt hole 
adjacent to the front flange. 


Cause 
The landing gear pivot shaft, 
Bendix P/N 66584, was 
broken approximately 4 in. 
from the strut. 


Reporting activity attributes 
failure to improper heat- 
treatment of steel in inner 
cylinder, resulting in crys- 
tallization. 


Overstress from continuous 
hard landings. 


Link assembly, main land 
gear upper retracting P/N 
5066434-1, broke at lower 
bolt hole. Believed due to 
previous hard landings. 

Indeterminate. 


Caused by a hard landing. 


Believe due to hard landings 
aggravated by sharp fillet 
radius at the base of the 
cylinder head. 


Indeterminate. The failed 
bolt, although forwarded by 
the reporting activity, has 
not been received by this 
Headquarters. 


Failure of fitting assembly, 
brace strut attaching, P/N 
2116632. 


Corrective Action 

This failure is believed caused by cracking 
the pivot shaft during a previous landing 
when either an extremely hard landing 
or a hard landing with full ammunition 
and fuel load was made. A design 
change is not believed warranted. 

Unsatisfactory reports of subject assembly 
have been negligible. Examination of 
defective parts, forwarded to this Head- 
quarters, revealed material and work- 
manship to be acceptable and that 
failures were due to overstressing result- 
ing from excessively hard landings and 
operation of airplanes on rough terrain, 
Further failures should be reported 
without delay. 

Replace defective part. All parts of 

« landing gear should be checked fre. 
quently and repaired or replaced as 
necessary. 

In view of the small number of similar 
failures reported and the fact that this 
type airplane no longer is in production, 
no further action is deemed warranted. 


First failure of this type to be reported; 
considered isolated; no action con- 
templated. 


Project No. MAC-3M-84. No action to 
strengthen or redesign the main landing 
gear struts is warranted in view of the 
small number of failures reported to date 
and since static tests have proved that 
the landing gear meets the design re- 
quirements in all respects. 

Project No. MAC-2F-42. Since there is no 
record of similar gear failures in the 
history of B-, no action relative to 
redesign will be taken. Replacement of 
the defective head cylinder with a similar 
serviceable part is the only action 
considered necessary. 

Project No. MAC-4C-1161M. Without 
examination of the bolt, the cause of 
failure cannot be ascertained. Previous 
hard landings may have caused failure, 
or seizing of the bolt due to inadequate 
lubrication may have been contributory. 

Project No. 16-1611-CM. Previous failures 
of fitting assembly, P/N 2116632, were 
attributed to lack of, or insufficient 
lubrication. Tolerances will permit a 
condition in which it is possible for the 
clearance between the fitting bushing 
and attaching bolt, P/N 1113280, to be 
as low as 0.0005 in. Freezing of the 
assembly and possible failure of the 
fitting could readily occur if the specified 
lubrication requirements are not com- 
plied with in periodic servicing of the 
assembly. 


Excessive loading due to a Project No. 16-1844-C. Examination of 


hard landing. 


(Table 2 continued on next page) 


the failed support landing-gear trunnion, 
P/N 117-33104, submitted with the 
report revealed no material weakness 
except for a slight deficiency in hard- 
ness. However, the margin of safety de- 
signed into the trunnion exceeded the 
loss of strength because of this dis- 
crepancy and should not have caused 
any trouble except when subjected to 
excessive loading. As the unsatisfactory 
condition reported is not believed to be 
of a general nature, replacement of the 
landing-gear trunnion support should be 
accomplished. 
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TABLE 2 (Continued) 


Date of 
U.R. Digest Model 
3/5/47 
Series 


Diffieulty 

Airplanes are being used for 
long-range weather recon- 
naissance and only make ap- 
proximately eight landings 
during 100 hours of operation. 
Reporting activity states that 
these circumstances should 
warrant extension of the 100- 
hour inspection of the landing 
gear and mechanism (Section 
10) to 300 hours. 


Failure of rear leg, landing gear 
retracting and attaching bolts 
part Nos. AN9-50 and 188011, 
allowing landing gear to col- 
lapse. 


4/20/47 C- 


because of vertical or side loads but rather because of 
loads in the drag direction either forward or aft. The 
most predominant source of failure appears to be in the 
drag bracing itself or its attachments or in some ele- 
ment of the structural system transmitting drag forces 
throughout the structure. Hard landings, on the other 
hand, are considered as such because of severe sinking 
speed, lack of flare, or excessive drifting. A case 
reporting the forward speed at contact, however, has 
yet to be found, and it is doubted if pilots or operators 
are aware that touch-down speeds have a profound 
effect upon landing-gear drag forces. 

Table 2 also includes some representative examples 
of landing-gear troubles which did not necessarily result 
in accidents. Such cases would seem to indicate that 
the landing-gear structure is experiencing loads in 
normal landings that are a large percentage of the 
landing gear ultimate strength. It is known from 
structural fatigue data that a rather large number of 
repetitions of loads just short of ultimate are necessary 
to produce damage or failure, and therefore those 
cases where failure is ascribed to prior damage must be 
viewed as cases where failure perhaps was imminent 
at the time initial damage occurred. 

Attention is again invited to the concise statement of 
the situation as reported by the A.A.F. Office of Flying 


Probably 
landings imposing strain in 
excess of the design load. 


Corrective Action 
Project No. 16-3313-AM. An investigation 
conducted by this Headquarters on 
landing-gear failures to date on B- air- 
planes revealed that considerable diffi- 
culty has been experienced with the 
component parts of the landing gear, 
resulting in extensive and major damages. 
The contributing factors to these numer- 
ous failures were insufficient or no 
breakaway torque, improper limit switch 
adjustments, malfunction of safety 
switches, possible misalignment of drag 
struts and drag links, defective tires and 
brakes, etc. In view of this, extending 
the 100-hour inspection to 300 hours is 
not considered advisable at this time. 
However, service tests are being con- 
ducted on long-range weather recon- 
naissance airplanes, and reports from 
these tests will be studied to determine 
the advisability of reducing the landing- 
gear inspection frequency on all B- 
airplanes and larger airplanes. 
excessively hard Project No. 17-554-B. Applicable Hand- 
books and Inspection and Maintenance 
Guides will be revised to direct the 
removal of the rear legs, landing gear 
retracting, at each 300-hour inspection 
period and the inspection for signs of 
wear of the legs and attaching bolts. 
Where the attaching bolts show signs of 
shear, the rear leg, as well as the bolt, 
should be condemned and replaced with 
serviceable parts. All rear legs, bolts, 
and bushings that show signs of wear 
should be replaced with serviceable parts. 
This inspection combined with the daily, 
25-hour and 100-hour inspections pres- 
ently specified should reduce the possi- 
bility of failure of these parts. 


Cause 


Safety (see p. 265). Sixty cases of landing-gear 
collapse in 90 days, of which 30 cases occurred during 
normal landing roll, means that once every 3 days in 
this period a landing gear collapsed in a ‘‘normal’’ 
landing. It should be noted also that this occurred 
only to aircraft within the continental United States 
and does not include accidents in combat theaters. 
These cases cannot be dismissed as landings where the 
gear should have failed. 

Certainly some of the many accidents resulting from 
landing-gear failure can be attributed to maladjust- 
ment, inadvertent retraction by the crew, etc., but far 
too many can only be ascribed to loads that are simply 
in excess of the ultimate strength and many of which, 
in turn, occur in normal landings. 


DESIGN CRITERIA AND CONDITIONS 


Landing gear and related structural design have 
been handicapped through the years by lack of a satis- 
factory means of expressing exactly what the landing 
gear should permit the airplane to do in landing and 
taxiing. Objective criteria have only recently been 
introduced in design specifications in the form of limit 
and ultimate sinking speeds and limit design touch- 
down speeds. Until these recent innovations, landing- 
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gear strengths have been prescribed in terms of arbi- 
trary load factors, arbitrary drag and side force coeffi- 
cients, and’ assumed airplane attitudes, and thus no 
criteria existed to determine what the airplane is sup- 
posed to be able to do without damaging the structure. 
This important element has been left up to the pilot 
and still is, to a great extent. 

In contrast to this, the strength of the structure for 
flight loads is completely specified by the familiar 
V-g diagram and by design gust intensities determined 
from atmospheric research. Furthermore, an airspeed 
indicator is provided to inform the pilot the magnitude 
of V and an accelerometer to indicate g. No such 
instrumentation is provided for the parameters im- 
portant in landing. Landing-gear load design philos- 
ophy continues to be to provide strength sufficient 
for what a pilot might do and yet expect the structure 
to survive. 

Consequently, landing-gear structural criteria consist 
of a list of “‘conditions” intended to include all reason- 
able circumstances, and the list continues to grow in 
length and severity. This increase in rigor of landing- 
gear structural requirements has undoubtedly sprung 
from the general awareness of landing-gear service 
failures, but the service statistics already shown do 
not reflect a corresponding improvement. Service 
records (see Tables 1 and 2) show that structural failures 
of the landing gear are predominantly in the drag 
direction and not necessarily from high vertical or side 
loads. Thus, drag loads for use in design have been 
increasing over the years (Fig. 5), but no attendant 
reduction in service failures from drag loads can be 
detected; in fact, Figs. 2, 3, and 4 indicate either an 
increase in gear failures during the past years or a 
trend at best unaffected by these more severe criteria. 

The design sinking speed for commercial and military 
(except carrier-based) aircraft has been about 12 ft. 
per sec. for some years. This means that the airplane 
is theoretically designed to sustain without failure 
an unflared landing of 720 ft. per min. Such a landing 
is so extreme that it is doubted that any pilot or 
observer would expect the airplane to sustain it without 
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failure. Yet the landing gear, and often the complete 
airplane, is drop-tested to this and even higher sinking 
speeds without failure and usually without any de. 
tectible damage. 

In Figs. 6 and 7 have been summarized some landing 
statistics contained in N.A.C.A. tests over an extended 
period with several different types of airplanes.’-” 
These tests were made to obtain, deliberately, landings 
covering a wide range of severity and, hence, must be 
considered to represent conditions that are certainly, 
if anything, “‘on the high side.’’ Out of a total of 244 
test landings, only two exceeded 11 ft. per sec., and 
these were both simulated blind landings. Out of 
173 flared landings, 115 were 3 ft. per sec. or less, and 
only 6 were 8 ft. per sec. or more. Even the unflared, 
blind landings averaged only about 5 ft. per sec., and 
only 8 out of 71 were 8 ft. per sec. or higher. 

The point to be made from the foregoing facts is 
that the design sinking speed of 12 ft. per sec. is indeed 
a “hard” landing, and it is suspected that many service 
failures ascribed to ‘“‘hard” landings are in reality con- 
siderably less severe than the 12 ft. per sec. unflared 
landing for which the airplane is designed and demon- 
strated to be adequate. 
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landings. 
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SpPEcIFIC ACCIDENT CASES 


The foregoing discussion has, of necessity, been 
general, since published service statistics contain little 
specific information about the circumstances surround- 
ing landing-gear accidents and still less regarding the 
diagnosis of the trouble or any action taken toward 
prevention of recurrence. The discussion does, how- 
ever, make clear the preponderance of landing-gear 
failures as a structural accident cause, the apparent 
paradox that many of such failures occur in normal 
landings and that design conditions are extremely 
severe and probably seldom realized even in so-called 
hard landings. In order better to illustrate the funda- 
mental problem with which this report is presently 
concerned, some specific cases of the Lockheed Con- 
stellation are discussed below. In these cases, rather 
complete information is available. 

It is conservatively estimated that, at this writing, 
Constellation airplanes have made approximately 
30,000 landings, of which a large proportion were in 
pilot training and familiarization usage. A good 
example concerns one air line that accumulated 1,366 
landings in 691 hours of flight time on a single airplane, 


"principally checking out pilots. Since Constellations 


are in use on several different air lines in both domestic 
and international operation, a wide variety of condi- 
tions and severity of landings have been experienced. 
There are many cases of landings where all individuals 
concerned have felt that the circumstances were so 
bad that complete gear failure could reasonably have 
occurred but did not. For example, there were these 
two: 

(1) A landing was made in which the wheels struck 
a sea wall in full flight. Contact was made several 
inches below the top of the wall as evidenced by tire 
tracks. A normal landing followed and subsequent 
inspection showed rather minor damage considering 
the circumstances. This damage consisted of partial 
compression failure of the drag strut, which is remark- 
able since aft drag loads produce tension in the drag- 
strut system. 

(2) A landing was made in which one wheel on 
one gear struck an obstruction (reported to be the 
stump of a tree) in full flight just prior to landing. 
This impact failed the torque arm and caused the gear 
to twist around 90° to the flight path. A normal land- 
ing was thereupon made in turf with the wheels on 
this side perpendicular to the path of the airplane. 
No damage resulted although these wheels plowed a 
deep furrow in the earth for a distance of about 
800 ft. 

These cases illustrate rather well the severity of 
landing which this gear is capable of sustaining, and 
doubtless similar experiences with other types of air- 
craft can be cited. In spite of these cases, however, 
there have occurred two instances of landing-gear 
failure in Constellation landings where the conditions 


AIRCRAFT LANDING GEARS 


(at least those controlled by the operator) were not 
particularly severe insofar as they could be determined. 
These cases were: 


(a) Air-line landing at Santa Maria Island (Azores), 
November 5, 1946. This landing appeared from all 
reports to be a moderately fast landing with moderate 
sinking speed, at a comparatively low gross weight and 
on a fairly abrasive runway that was wet at the time. 
After a run of about 1,500 ft., the left main gear 
collapsed rearward. No personnel injuries or fire 
resulted, but considerable material damage occurred. 
The airplane was subsequently repaired and returned to 
service. Complete investigation of this accident 
established that initially a compression failure of the 
drag strut occurred and that this was ultimately 
followed by complete tension failure of the damaged 
strut during brake application. 

(b) Air-line landing at Logan Field, Boston, May 
3, 1947. This landing was made by particularly 
experienced crew and appears from all reports to have 
been moderately fast with moderate descent rate, at 
normal gross weight and on a brand new smooth con- 
crete runway somewhat wet at the time of landing. 
After a run of about 2,000 ft., the right main gear col- 
lapsed rearward. No personnel injuries or fire resulted, 
but considerable material damage occurred. This air- 
plane has subsequently been repaired and restored to 
service. This landing caused the complete failure of the 
wing bulkhead structure (to which the drag strut at- 
taches) under the action of a fension load in the drag 
strut. No important failure of the gear itself occurred. 
Initial failure was determined by all investigators to 
have occurred in a particular lug in the drag bulkhead 
which is the precise location of first failure in destruc- 
tion static tests conducted by the A.A.F. at Wright 
Field at ultimate load. 


In the above two cases, the failures were duplicated 
in static tests and occurred at or above ultimate loads, 
which, as briefly discussed earlier in this report, are ex- 
tremely high. (Further discussion of landing-gear 
strengths is found later in this report.) No evidence of 
defective materials or workmanship could be found. 
In both cases, visibility, ceiling, and wind conditions 
were satisfactory or favorable. 


Thus, here were two landings where the conditions 
appeared not unduly severe and most probably did not 
even approach the conditions designed for, and yet 
loads in the structure exceeded ultimate values. Con- 
sequently, there appears to be some condition or phe- 
nomenon present under certain circumstances which 
causes excessive structural loads, and this phenomenon 
must be found in order to explain the cases cited above 
and also those cases scattered throughout the service 
records of commercial and military aircraft where struc- 
tural failures occurred in normal landings and were 
ascribed to some inconsequential cause or were unex- 
plained entirely. 
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POSSIBLE EXPLANATIONS OF FAILURES 


In the following sections, several items that would 
offer conceivable explanations for landing-gear failures 
in normal landings are discussed. These explanations 
haveJbeen advanced from time to time, and each has 
considerable grounds for suspicion as the source of gen- 
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eral landing-gear structural trouble. It is not possible 
in this paper to discuss completely each item, but suf. 
ficient facts are at hand to establish the relative mag- 
nitude and importance of each theory. 


(1) Landing Gears Inherently Too Weak 


This is a fairly obvious explanation, and it would 
seem that there should be few failures that could not be 
prevented by heavier and stronger members, fittings, 
and attachments. However, as previously pointed out, 
design loads have been increasing over the years, and 
landing gears have, as a result, become heavier and 
stronger but with no attendant reduction in service 
failures. Service records also show that landing-gear 
failures do not predominantly involve types designed 
for less severe criteria. Fighters are, in general, de- 
signed for landing loads about twice those for other 
types, and yet their record is no better. It has been 
British design practice for some time to use an ultimate 
landing load factor of 4.00 for all types including fight- 
ers, while this factor has been 7.00 for American fight- 
ers. British design to the lower strength factor would 
certainly not have continued if their record of gear 
structural failures was appreciably worse than that for 
American fighter airplanes. 

The Constellation landing gear, for example, is ex- 

tremely strong and meets C.A.A. requirements with 
ample margins, as do other commercial aircraft. The 
strength of this gear may be compared with that of 
other aircraft in a simple manner by means of a dia- 
gram known as a “Landing Gear Strength Envelope.” 
This diagram is described in detail in reference 18 and 
briefly consists of a plot of the vectors of the actual 
forces required to produce failure in the structure using 
the axle in the fully extended position as the origin. The 
locus of the extremes of these vectors forms an envelope 
or boundary enclosing all possible combinations of ulti- 
mate forces in a particular reference plane. Fig. 8 is 
such a diagram for the Constellation landing gear, show- 
ing actual strength as determined in various destruc- 
tion tests. Comparison is made in Fig. 9 where actual 
landing-gear strength divided by landing weight is 
plotted for several different airplanes. This shows that 
the Constellation compares favorably in strength to 
these other airplanes—-commercial transports and 
bombers. The effect of service change (S.B. No. 287, 
Increasing Bolt Size in Upper Drag Strut) is illustrated. 
The strength of the gear for forward loads was increased 
from a value comparing well with these types to a new 
value about 35 per cent higher than before. This 
change was rather recent, however, and the bulk of the 
landings amassed to date were on landing gears having 
the former strength. 

Increased strength also must be provided generally 
throughout a structure or else troubles are merely moved 
to new locations. The structure resisting drag loads is 
analogous to a chain where strengthening the weakest 
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AIRCRAFT LANDING 


link raises the general strength level only to that of the 
next weakest element. This point is not a matter of con- 
jecture, for experience on PV-1 aircraft bears this out as 
does an examination of service reports on other aircraft. 

Furthermore, it would appear that the existence of 
ultimate loads in ‘“‘normal’’ landings requires increasing 
the general strength level for landing loads by a large 
factor, perhaps the order of two. This is an enormous 
weight penalty and should raise the weight of the land- 
ing gear and affected structure to about 15 to 20 per cent 
of the empty weight, which is intolerable. 

Consequently, it appears that increased landing-gear 
strength as a solution to recurrent landing-gear failures 
would at best be doubtful and would in any event cost a 
considerable penalty in weight. 


(2) Fatigue Effects 


There is considerable interest in fatigue effects and 
the inherent service life available from a given structure. 
Fatigue failures are definitely in evidence in trouble re- 


- ports, although much has been accomplished through 


improved detail design and elimination of stress con- 
centrations. However, fatigue failures have occurred 
in many instances after such short service periods that 
they can only be due to a comparatively small number 
of repetitions of extremely high loads. For a fatigue 
crack to occur after, say, 5,000 hours of operation means 
that perhaps 1,000 landings have precipitated failure. 
Fatigue literature shows that practically any structure 
will sustain without damage several thousand repeti- 


_tions of loads that are two-thirds to three-quarters of 


the ultimate static strength. Thus, loads considerably 
in excess of limit design values and probably closely ap- 
proaching ultimate must be experienced for fatigue 
cracks to develop in short service periods. The two 
cases of Constellation gear failure revealed no evidence 
whatever of fatigue damage. Other accident records 
show, also, too many cases where no evidence of prior 
damage existed at all. This is probably to be expected 
in view of the rigorous and frequent overhaul, inspec- 
tion, and replacement that landing gears receive, with 
the result that cracked parts are caught before serious 
trouble can occur. 

Thus, while fatigue can be cited as contributing to 
the poor service record of landing gears, it must be 
considered an effect rather than a cause, and, in fact, the 
very existence of fatigue cracking after short service 
points to the existence of abnormally high loads in rou- 
tine landings. 


(3) Shock-Strut Behavior 


Most landing gears in recent years consist of an ap- 
proximately vertical shock strut to resist vertical loads 
and braced for drag and side loads in a manner that 
subjects the shock strut to bending. The piston of the 
shock absorber is invariably the lower member, which 
moves into the cylinder located above and rigidly sup- 
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ported. The shock absorber itself is usually of the air- 
oil type, depending on air under pressure as a taxiing 
spring (augmenting the tire) and to extend the piston 
and depending on the oil for absorption of large energy 
input and rebound damping. 


There have appeared several difficulties with this 
arrangement which have been cited as reasons for un- 
necessarily high loads in normal landings. One of these 
is that design for the best compromise between stroke 
and loads results in a demand for high strut efficiency, 
frequently leading to a metering pin to maintain the 
load in the “‘design’’ landing condition within accept- 
able limits. This design results in utilization of less than 
full stroke in landings less severe than the limit condi- 
tion, the reduced energy being absorbed by higher loads 
in combination with lesser piston movement. This 
characteristic is inherent in any conventional air-oil 
shock strut, with or without a metering pin, and is il- 
lustrated in Fig. 10. A most desirable characteristic, 
which the author mentions as a challenge to designers, 
would be a shock absorber that utilizes the full stroke 
available regardless of the landing severity, thereby 
reducing the loads in light landings to values more 
nearly proportional to the quality of the landing. 


Furthermore, air pressures determined by required 
piston taxi positions are high when the piston is fully 
extended (also illustrated in Fig. 10). These pressures 
are established by permissible bottoming forces and by 
forces required to cause the wheel to follow irregulari- 
ties in ground contour. However, these fully extended 
air pressures must be overcome before movement of the 
piston occurs, and thus vertical loads in ‘‘soft’’ landings 
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are increased by virtue of this characteristic. Bending 
of the shock strut itself also occurs with attendant cyl- 
inder ovalization and, together with friction forces due 
to side and drag loads, necessitates additional vertical 
loads to produce movement of the piston. This latter 
condition is, of course, made worse if drag or side loads 
are abnormally high; and the smaller the shock strut 
travel, the greater is the moment arm of these forces to 
the supporting structure. 

The above features are rather common complaints 
about the modern aircraft shock absorber, and there is 
much room for improvement which can only come from 
research and experiment. This observation is corrobor- 
ated in some German work.*! However, shock absorb- 
ers and frequently complete airplanes are subjected to 
drop tests in which design vertical and drag loads are 
simulated and excessive loads in the design conditions 
are not tolerated. Less severe conditions result in loads 
that, although not proportionately lower than those in 
extreme conditions, are certainly not higher. Thus, 
these inherent strut characteristics cannot be viewed 
as a primary explanation of destructive loads in a nor- 
mal landing, although such features undoubtedly con- 
tribute to the generally high level of loads that induce 
early fatigue failure. Of course, as mentioned above, ex- 
cessive drag (or side) loads can amplify these charac- 
teristics, and therefore some further discussion of this 
point is presented in subsequent paragraphs. 


(4) Excessive Ground Friction 


Design drag coefficients have been increasing during 
the past years, as was illustrated in Fig. 5. For many 
years, the design aft drag load was 25 per cent of the 
design vertical load. Considered inadequate, presum- 
ably because of service experience, this figure was raised 
in 1941 to 0.33, an increase in drag strength of 33 per 
cent. Friction coefficients for ‘braked landings’’ and 
brake torque requirements have been consistently 
around 0.55. These coefficients have been used in gen- 
eral aircraft design. However, automotive literature for 
some time has shown that friction coefficients between 
tires and pavement may have much higher values, 
reaching 0.80 and even greater than 1.00. An effect 
with speed is noted, and, for slow speeds between 
smooth rubber and such surfaces as glass and steel, co- 
efficients as high as 4 and 5 (400 and 500 per cent of the 
normal force) have been measured at the Bureau of 
Standards.* These latter values are perhaps of aca- 
demic interest only, but they illustrate the range of 
values which a friction coefficient may exhibit. 

Consequently, if ground friction coefficients of around 
0.80 can be realized, and aircraft landing gears have 
been designed until recently for only 0.25 to 0.33 times 
the vertical, it would at first appear that here is a source 
of loads in excess of design. However, it should be 
realized that such coefficients are significant (insofar as 

drag loads are concerned) in two conditions—viz., dur- 
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ing brake application and during the wheel spin-up in. 
terval. In the former case, vertical loads are little 
higher than 1g (design limit = 1.33g) and the shock 
strut is in a favorable position (static). In the latter, 
the wheel skids only until it reaches ground speed, and 
the maximum value of the ground friction force is there. 
fore limited by the value of the vertical force (multiplied 
by friction coefficient) at this instant. Generally, 
wheels come up to speed in a fraction of the time re. 
quired to reach peak vertical load, and thus the ground 
friction force is usually much lower than the maximum 
vertical force. Now an inspection of the strength dia- 
grams in Fig. 9 shows that these aircraft can all sustain 
drag ratios of around 0.8 to 1.0 up to substantial verti- 
cal loads and, because the shock strut compresses 
under vertical loads, can sustain still higher drag 
forces. (Diagrams are based on fully extended struts 
in all cases.) Furthermore, some of the cases of failure, 
at least, have occurred under conditions where the 
friction is known to be lower than in other uneventful 
landings. For example, the two Constellation inci- 
dents cited earlier both occurred on wet (but not 
flooded) runways where friction coefficients are known 
to be lower. Similarly, Army and C.A.A. records show 
accidents of this type on runways known to have com- 
paratively low friction characteristics. 

Thus, high coefficients of ground friction cannot be 
considered a satisfactory explanation of landing-gear 
failures in ostensibly normal landings but rather must 
be viewed as another aggravating and contributing 
element to generally high landing forces and as a land- 
ing variable that has a wide range of values. 


(5) Structural Deflections 
As the landing gear and its supporting structure de- 


‘flects under the action of the external forces, secondary 


effects and geometry changes result. Such effects may 
increase the loads in structural elements over those cal- 
culated for the undisturbed geometry and may become 
significant in certain gear geometries and flexible 
structures. In general, however, these effects are rather 
small, seldom amounting to 10 per cent, and in any 
event are provided for either in the stress analysis or in 
structural tests. The airplanes whose gear strength 
diagrams are shown in Fig. 9 have all been tested to ul- 
timate loads in full-scale tests, and most of those cited 
in the accident cases of Tables 1 and 2 have been simi- 
larly substantiated. Consequently, secondary effects 
due to structural deflections do not appear to be a ma- 
jor consideration, although they cannot be ignored as a 
source of structural loads in excess of nominal values. 


DISCUSSION OF DEDUCED CAUSE 


Previous discussion has shown, it is felt, that while 
there may be several aggravating or contributing causes, 
there has yet to be found a condition of sufficient con- 
sequence to explain why aircraft can make thousands of 
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satisfactory landings under the most adverse conditions 
and still experience landing-gear failures in landings 
that appear to be not unduly severe. An underlying 
cause, however, of sufficient magnitude to explain a 
great many structural gear failures of the past has been 
found. This cause has been deduced from scattered 
evidence from many sources, some of which has been 
presented in this paper. In subsequent paragraphs, 
this phenomenon is discussed in some detail, together 
with its applicability to the cases available for study. 

The landing gear is an elastic structure and has a 
natural frequency of vibration in a mode corresponding 
to fore-and-aft motion of the axle. This oscillatory 
motion has for some time been known to occur under 
the action of wheel spin-up drag forces that accelerate 
the wheels up to ground speed. Some of the first quan- 
titative work on this subject is reported by Bergen*? 
and Howland.** There have also appeared from time to 
time the results of tests on other airplanes by the N.A.- 
C.A. and the A.A.F. (for example, references 19 and 34). 
These papers show that it is a fundamental character- 
istic of these modern airplanes to experience fore-and- 
aft oscillations of the gear as a result of wheel spin-up 
forces. Fig. 11 is a typical plot of forces measured in an 
actual landing on a four-engined bomber!’ which illus- 
trates this vibration and some of its effects. 

Lockheed has conducted tests in which actual meas- 
ured loads were obtained in landings of instrumented 
airplanes, including the Constellation, and has ob- 
served the same phenomena. The results of some of 
these measurements on the Constellation are shown in 
Fig. 12 where the measured drag force vs. time curves 
obtained in several landings have been superimposed 
with their common point at0. The traces to the right of 
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0 have been adjusted to the average vibration frequency 
for clarity, but the traces to the left are exactly as 
measured. 

Now the wheel spin-up drag force which excites this 
vibration is a function of many variables, some of which 
have been touched upon in previous discussion. Some 
of these variables are: (1) sinking speed and flare; 
(2) forward speed and wind; (8) attitude of airplane; 
(4) tire inflation; (5) tire tread and its condition; 
(6) runway—material, condition, and deflection; (7) 
variation of friction coefficient with speed, etc.; (8) 
shock strut behavior (vertical force and rate, and ex- 
tension); (9) gross weight of airplane; (10) brake 
drag and other variations in wheel inertia; and there 
are undoubtedly others. In addition, many of these 
items are interrelated. For example, vertical load ef- 
fects the tire rolling radius, which in turn effects the 
r.p.m. that the wheels must achieve to match the ground 
speed. These variables not only affect the magnitude 
of the wheel spin-up ‘drag force but also its duration— 
its rate of application. 

It will be noted that of these many variables only the 
first three are under control of the pilot, and these are 
the items that determine whether or not a particular 
landing is good or bad. Furthermore, all these items 
have a wide range of values. The friction coefficient 
can vary from around 0.02 for an icy field to perhaps 0.8 
or 1.0 for dry concrete or smooth macadam. The gross 
weight at landing can easily vary by 25 per cent among 
long-range, short-range, and training loads. Brakes 
have a wide scatter in free-wheeling torque because of 
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adjustment, wear, dirt, etc., and shock-strut sensitivity 
to several conditions has already been discussed. Thus, 
it appears that, while the parameters controlled by the 
pilot may have nominal values, there are a great many 
other factors present over which he has no control and 
which have a wide range of values. 

Forces applied rapidly to an elastic structure produce 
deflections and consequent stresses quite different from 
those resulting from the same forces applied in a static 
or slow manner. The resulting deflections are higher 
than the static movements by an amount depending 
on the relationship between the rate of application of 
the force and the natural frequency or response of the 
structure. This is discussed in many vibration texts 
and is treated clearly in a recent paper ‘‘On the Effect 
of Sudden Variation in Force” by Herrey.?? The phe- 
nomenon is usually referred to as “‘dynamic overloading” 
or transient stresses, and Fig. 13 is a reproduction from 
that paper to illustrate this dynamic response of an elas- 
tic structure to a rapidly applied force. The reference 
paper also brings out simply that the phase position 
of the deflected structure at the instant of peak force 
has an important effect and that the magnification 
factor caused by the dynamic response can reach a lim- 
iting value of two. These dynamic forces are the forces 
that have been measured in actual landings, and thus 
ground forces (or forces at the axle) deduced from them 
are “apparent” and not actual. Thus, apparent ground 
friction coefficients in excess of 1.00 are most likely the 
result of this phenomenon. 

Data on natural frequencies of landing gears are not 
generally available, but measurements taken on a few 
airplanes radically different in size, geometry, and tac- 
tical use show that these natural frequencies are around 
6 to 20 cycles per sec. or natural periods of about 0.16 to 
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0.05 sec. Measurements of wheel spin-up forces indj- 
cate spin-up times that lie in this range; therefore, 
critical dynamic response is possible. 

Consequently, since landing gears have a natural mode 
of vibration of the same order as wheel spin-up times and 
since wheel spin-up rates are functions of many uncon- 
trolled variables each having a wide range of values, there 
then exists a random chance that these values may occur in 
combination such that dynamic overloads of large magni- 
tude can be experienced on any landing gear. Further- 
more, a condition may exist where some other portion 
of the structure responds to the landing-gear forcing 
vibration and thus experience dynamic repetitive 
loads. 


Here, then, is a source of excessive loads in normal 
landings of sufficient magnitude to explain a consider- 
able number of failures. Since the dynamic overload- 
ing can have different values, various kinds of failure 
can occur. Lacking sufficient data on other aircraft 
the Constellation experience will again be used as the 
example. Since the landing gear is (or was at the time) 
stronger in the aft direction than the forward, a dynamic 
overload can be experienced which is not sufficient 
to cause tension failure at the initial peak but which is 
sufficiently high to cause compression failure at the for- 
ward rebound. This is believed to complete the expla- 
nation of the accident (Case a) cited above. Likewise, a 
higher overload can be conceived such as to cause imme- 
diate tension failure or important damage at the initial 
rearward peak load which reduces the strength of the 
drag system in tension to the point where subsequent 
braking forces destroy the structure. This is believed 
the true explanation of the accident (Case b) reported 
above. (It is possible that more rapid rates of drag-load 
application occur on wet surfaces such as existed in 
both of these cases. Under light vertical loads, it ap- 
pears that the tire may skid on a film of water with low 
friction, suddenly changing to a higher value as the 
increasing vertical force causes the tire to ‘‘squeegee” 
the water from beneath it.) 


When excessive loads of this nature are experienced, 
damage to the structure radically changes the frequency 
and produces considerable damping. Hence, immedi- 
ate, complete failure would in most cases not occur. 
Subsequent use of brakes to stop the landing roll, how- 
ever, entails substantial loads upon the drag system, 
which loads are in no way reduced by prior damage, 
and thus the complete failure so often observed at this 
point in the landing operation would logically occur. 


To this phenomenon also can now be attributed the 
relatively high loads necessary to cause many of the 
fatigue failures after short service which appear so 
frequently in the service records of commercial and 
military aircraft of all types. It is suspected also that 
other troubles elsewhere in aircraft structures perhaps 
may be ascribed to the dynamic response of these struc- 
tures to the heavy landing-gear forcing vibration. 
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DISCUSSION OF OTHER SOURCES OF VIBRATORY MOTION 


Although there are no cases known to the author, 
there must be admitted the possibility of encountering 
irregularities in runway surface of the proper spacing 
to cause resonant response of a landing gear or a motion 
analogous to aerodynamic flutter. This condition 
would appear to have a low probability, since precise 
“tuning’’ of speed, natural frequency, and bump spac- 
ing would be required to excite it. ' 

There is, however, a phenomenon peculiar to wet 
runways that merits a more detailed discussion. In 
Fig. 14, there are presented” curves of skidding friction 
coefficients versus speed of slip for automobile tires on 
various pavements. Reference 26 contains similar 
data. It will be noted from these curves that the coef- 
ficients of friction for dry surfaces are high and rela- 
tively constant with skidding speed, while those for wet 
surfaces decrease markedly as the skidding speed in- 
creases. Other data (not shown) for ice, sleet, and 
packed-snow-covered surfaces show low coefficients 
that remain constant with speed of slip. 

Now, if the wheels are skidding (locked) on a wet 
surface, the condition illustrated in Fig. 15 can readily 
be conceived. Under the steady friction force, the ini- 
tial conditions of gear deflection and skidding velocity 
(Vo) are as shown by the left-hand portion of these 
curves. Upon encountering a bump, a dry spot, a 
more abrasive spot, or a puddle on the runway, an in- 
creased drag (friction) load results, the gear deflects fur- 
ther rearward, and the relative skidding velocity is de- 
creased to V;. When the disturbance is passed, the 
friction force decreases, the gear deflects forward, and 
the skidding velocity increases to Vo. Thus it will be 
seen that the change in friction with skidding speed 
which occurs on wet surfaces results in an oscillatory 
ground force capable of producing resonant vibration, 
since the friction force is always in phase with the velo- 
city of the gear motion. It would appear that these 
vibrations would be divergent until the relative skid- 
ding velocity reached zero, until a region of stable fric- 
tion was reached, or unless adequate damping existed. 
There is some evidence that a few of the failures re- 
ported in the past have resulted from this phenomenon. 


DISCUSSION OF Nose LANDING-GEAR STRUCTURES 


In the case of nose landing gears, significant dynamic 
response at spin-up can exist but with much lower prob- 
ability than that for main gears because the variables 
involved are fewer and have narrower ranges. In the 
vast majority of landings, the nose wheel is not touched 
until the forward speed has reduced, and, when it does 
so, it can make contact at substantially one attitude 
only. Generally, then, the nose gear is subjected to 
vertical loads that result from pitching motion only, 
since the main gear has already absorbed all vertical 
velocity. Not having brakes, the nose wheels are not 
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subject to variable brake drag and, being castored, 
cannot develop significant side loads in normal landings. 
Furthermore, since the pilot is usually located almost 
directly over the nose gear, he has a much more accur- 
ate ‘‘feel’’ of the forces and distances there; hence, he 
can and does generally succeed in subjecting the nose 
gear to comparatively light forces in landing. The ab- 
sence of brakes also eliminates the possibility of the 
skidding oscillations discussed above. 

In addition, it should be pointed out that the most 
critical design conditions for the nose gear often occur 
in ground handling, such as the use of symmetrical 
and unsymmetrical power against a blocked nose wheel, 
turning, and stopping conditions. These conditions 
usually cause the nose gear to have considerably more 
strength than would be dictated by landing conditions 
alone. Thus, it is concluded that the probability of 
experiencing damaging dynamic loads in the nose gear 
is exceedingly remote, and this is borne out by excellent 
service experience on those nose gear structures of which 
the author is cognizant. 


DISCUSSION OF POSSIBLE SOLUTIONS 


There are a few design recourses that offer possible 
solutions to the dynamic response of landing gears to 
the sudden application of wheel spin-up forces, and 
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these will be discussed briefly. These possible solutions 
are: 


(1) Increased Strength 


Increased strength has already been discussed earlier 
and found to offer at best a solution of doubtful success 
and expensive both in weight and cost. Aircraft land- 
ing gears are already extremely strong and expensive, 
and landing forces already design much of the aircraft 
structure. Furthermore, it will be seen from the fore- 
going discussion that greater strength cannot eliminate 
the probability of excessive dynamic loads unless the 
structure is made strong enough for approximately twice 
the most severe spin-up load possible.. This is tanta- 
mount roughly to doubling the weight of landing gears 
and supporting structures. Besides, oscillatory motion 
of the landing gear will still occur with the attendant 
possibility for excitation of other portions of the air- 
plane, and resonant response to runway irregularities 
and other vibrations is still possible. 


(2) Increased or Decreased Natural Frequency 


A raising or lowering of the natural frequency of the 
landing gear could improve or make worse a given in- 
stallation. The range of wheel spin-up forces and times 
is so great that it would appear that such change would 
merely result in some other critical spin-up combina- 
tion. Raising the natural frequency would appear the 
more likely of success because then only the most rapid 
spin-up loads would excite critical response. However, 
the natural frequency increases as the square root of 
the stiffness, and thus to double the natural frequency 
requires four times the stiffness. Furthermore, the 
frequency decreases as the square root of the mass, 
and, hence, a uniform strengthening of a structure with 
its attendant weight results in a disproportionately 
small chang? in its natural frequency. A radical change 
in geometry of landing-gear structures would offer the 
most beneficial change, but here again the frequency 
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can only be raised to that of the structure that supports 
it. The important mode of the supporting structure jg 
usually wing torsion, and this frequency is usually not 
much different from that of the landing gear. Further. 
more, regardless of the final frequency achieved, there 
remains the possibility of exciting important vibrations 
in other portions of the airplane and resonant response 
of the gear to external forcing vibrations. 


(3) Wheel Prerotation 


Wheel prerotation is one of the most promising means 
of eliminating the source of spin-up dynamic loads, and, 
although one of the “most-invented” items in aircraft 
design, it has reached only primitive development, 
There are many difficult problems connected with pre- 
rotation, and it is believed that, at this writing, only a 
few experimental ‘“‘prerotated’’ landings have been 
made, at least in the United States. However, mechani- 
cal and electrical means are known, and continued 
development will undoubtedly result in satisfactory 
prerotation. On the other hand, there are problems of 
reliability and secondary matters generated by prero- 
tation that merit discussion. 

From simply theory, the wheel spin-up drag force 
increases with the square root of the touch-down speed, 
other things being equal. (For example, see reference 
30.) This is graphically represented in Fig. 16, where 
it will be seen that 75 per cent prerotation is necessary to 
effect a 50 per cent reduction in drag force. This is 
based on a constant friction coefficient and the assump- 
tion that all other variables are likewise constant. When 
the change in friction coefficient with speed is accounted 
for, an even higher per cent prerotation may be neces- 
sary. 

British work on automobile tire friction®® shows 
the friction load between tire and pavement to increase 
rapidly from zero slip (100 per cent prerotation) to its 
full value at 20 per cent slip (80 per cent prerotation). 
Thus, a wheel prerotated to something less than 100 
per cent of ground speed may experience loads very 
nearly as high as when unprerotated, and the reduced 
momentum of the spinning wheels may simply result in 
shorter spin-up times. It should be noted also that 
20 per cent slip is possible in the case of ‘‘100 per cent” 
prerotation because of the effect of wind and of the 
variable rolling radius of the tire. No experience is at 
hand on the effect of wheel unbalance during the pre- 
rotation phase of the landing approach, but it’ would 
appear that more rigorous control of wheel balance 
would be necessary to avoid objectionable vibrations 
during this interval. Prerotation, of course, would have 
no effect upon other sources of vibratory motion. 
Therefore, although wheel prerotation appears to be a 
real solution to a major part of this fundamental prob- 
lem, it is not without many attendant difficulties that 
must be surmounted by future development and re- 
search. 
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(4) Damping of Fore-and-Aft Motion 


Serious efforts to damp the horizontal motion of the 
gear have not been found in the search made in this 
study. Reference 35 contains a brief description of a 
rubber device used for the drag strut of an Army two- 
engined trainer which experienced recurrent drag sys- 
tem failures during the war and of some landing meas- 
urements made therewith. Results indicate some ben- 
efit, and subsequent service was apparently improved. 
The degree of improvement, however, could only have 
resulted from the change in natural frequency, since 
rubber provides damping only by virtue of its hystere- 
sis, and its response is little, if any, different under 
widely different deflection velocities. Nevertheless, the 
device appears to have done some good and succeeded 
in reducing by about 40 per cent the number of oscilla- 
tions per landing, although peak rearward loads ap- 
pear to have been generally increased, which might be 
expected. 

The German literature indicates that some considera- 
ation was being given to this matter prior to V-E day. 
To date, no definite solutions have been found, but 
reference 31 indicates some experimentation with the 
shock strut to improve this characteristic on the ME 
262 on which “‘. . . breaking of a shock strut occurred at 
V,, = 200 km/h...” or a touch-down speed of 130 
m.p.h. 

If the vibratory motion of the landing gear could be 
hydraulically damped, however, dynamic overloading 
could be substantially eliminated, forward rebound 
loads reduced to comparatively smail proportions, and 
the number of oscillations per landing reduced from the 
ten to 20 severe reversals usually experienced on pres- 
ent-day aircraft down to one moderate cycle. This 
would eliminate the objections to the foregoing solutions 
and offer an enormous benefit to the fatigue life of the 
structure besides. Essentially, it is desired to provide 
damped motion of the axle in the drag direction (both 
rearward and forward) and to do so without reducing, 
and preferably improving, the effectiveness of the shock 
absorber provided for the vertical motion. This would 
then mean damping “‘in a plane’’ rather than damping 
“along a line,” as provided by present arrangements 
including the “‘linkage’’-type installation. There is thus 
required a dashpot to provide damping for rearward and 
forward motion of the axle of sufficient stroke and char- 
acteristics to cause fewer, and preferably no, subsequent 
vibrations. 


A shock strut located in the drag system would 
appear to offer this solution, although geometry and 
permissible travels might prevent optimum results from 
being achieved in some cases. A device similar to 
the shock strut seems suitable, for, although it has 
several shortcomings as discussed earlier in this paper, 
it is mechanically reliable and, when properly installed, 
damps the vertical vibration and bounce that would 
most certainly occur without it. 
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A device that damps fore and aft oscillatory forces 
has the further advantage over the previous solutions 
discussed of suppressing landing gear vibrations from 
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other sources besides the wheel spin-up condition. Vi- 
bratory drag forces due to brake chatter or to skidding 
chatter, while perhaps not possible to eliminate by a 
drag strut damper, should be prevented from build- 
ing up and their effect upon the structure significantly 
reduced. 


Test RESULTS 


As a result of this study, a hydraulic damper was de- 
signed for the drag system of the Constellation landing 
gear and subjected to numerous tests to determine its 
effectiveness.*® A diagram of this device is shown in 
Fig. 17. Essentially, the damper is a double-acting pis- 
ton with a mechanical spring for initial positioning. A 
metering pin is provided to control the motion for opti- 
mum energy absorption. The installation of the drag 
strut damper is illustrated in Fig. 18. The damper 
stroke permits a rearward movement of the wheels of 
about 6 in. from the initial forward position. 

The complete landing gear was mounted in the drop- 
testing machine, shown in Fig. 1, and was subjected to 
a systematic series of tests with the damper installed 
and with a rigid drag link. In these tests, ground forces, 
as well as structural forces, were measured by means of 
strain gages and a recording oscillograph, controlled 


sinking speeds were achieved through simple simulation 
of wing lift, and the wheels were rotated in reverse sense 
so that stopping forces at contact duplicated normal 
spin-up forces. 

Figs. 19 and 20 show some typical results of these 
tests. Fig. 19 is a plot of the load in the drag brace and 
the axle horizontal deflection as a function of time for 
the condition of 10 ft. per sec. sinking speed and a for- 
ward speed of about 95 m.p.h. In this case, the initial 
tension peak was reduced 32 per cent, and the compres- 
sion rebound was reduced to less than one-quarter its 
original value through the use of the damper; subse- 
quent oscillations: were reduced to negligible propor- 
tions. This case corresponds closely to the ‘‘limit’’ land- 
ing condition. Fig. 20 shows similar data for a more 
likely, though severe, condition of 95 m.p.h. forward 
speed and 6 ft. per sec. sinking speed. In this case, the 
initial tension peak was reduced to about one-third its 
former value, rebound was eliminated, and the motion 
was clearly dead-beat. In these and similar cases 
tested, comparable results were achieved. Of consid- 
erable significance is the fact that the damper reduced 
the oscillations from several severe reversals of load to 
one moderate half-cycle, which should unquestionably 
result in a major improvement in the fatigue life of the 
structure affected. Flight tests conducted on an air- 
plane instrumented with strain gages corroborated these 
results, and a marked reduction in the drag shock or 
“lurch” was observed by the occupants. 

A series of tests was also made on an instrumented 
airplane in an attempt to produce the phenomenon 
(previously discussed in this paper) of skidding oscilla- 
tions on a wet runway and to determine the effect of the 
drag strut damper upon them. Approximately 1,000 
ft. of runway surface (macadam) was wetted with fire 
hoses, and the airplane was skidded at various speeds 
throughout this length. The oscillations and tendency 
for unstable motion were indeed observed. Fig. 21 isa 
plot of one of these tests showing drag loads as a func- 
tion of time under similar conditions with and without a 
drag strut damper. The effect of the damper is clearly 
evident, reducing the magnitude of the loads enor- 
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mously and eliminating any divergent tendency. It 
will also be noted from this figure that a small fraction of 
the damper stroke was required to control the motion. 

It is beyond the scope of this paper to discuss in de- 
tail all the tests conducted and the results obtained in 
connection with the drag damping work carried on at 
Lockheed. The above results are intended merely to 
illustrate that hydraulic damping of horizontal gear 
motion appears to be a most effective means of eliminat- 
ing those sources of excessive structural loads that have 
caused a principal share of structural failures through- 
out the aircraft industry. 


CONCLUSIONS 


The conclusions reached in this study which have a 
bearing on the matter of landing-gear structural failure 
are these: 

(1) The landing gear of commercial and military air- 
craft has been responsible for the majority of structural 
failures and is second only to the power plant as a source 
of accidents and troubles insofar as “‘material’’-type 
failures are concerned. 

(2) Such structural failures are predominantly from 
drag loads in spite of steady increases in design drag 
loads. 

(3) A large number of these structural failures occur 
in ‘‘normal’’ landings. 

(4) Structural design conditions appear to be amply 
severe and are perhaps seldom realized even in ‘“‘hard” 
landings. 

(5) Landing-gear strengths of comparable aircraft 
are approximately equivalent and are all high. 

(6) Aircraft landing gears all appear to exhibit the 
same general behavior characteristics. 

(7) The predominant cause of landing-gear struc- 
tural failure is that loads are in excess of design values, 
even in “normal’’ landings, and abnormally high loads 
are frequently repeated in service. 

(8) The general strength level, fatigue, shock strut 
characteristics, extreme ground friction, and structural 
deflections aggravate and contribute to the problem but 
fail to explain service failures. 

(9) The fundamental cause is deduced to be dynamic 
overloading of landing gears resulting from critical 
rates of application of wheel spin-up drag forces. 

(10) Critical rates of application of wheel spin-up 
drag forces have a random probability of occurrence on 
any aircraft because of the number and wide range of 
uncontrolled variables that determine these rates. 

(11) Resonant gear response to forcing vibrations 
from the ground appears possible, and skidding on wet 
runway surfaces appears to be one likely source of such 
motion. 

(12) Nose landing gears have a lesser probability of 
such dynamic overloading because of inherently greater 
strength, of reduced number and range of variables 


affecting the spin-up forces and rates, and of the ab- 
sence of brakes. 

(13) Adjustment of landing-gear natural frequencies 
and wheel prerotation offer possibility of improving 
this condition but may not eliminate it entirely. 

(14) Hydraulic damping of the drag motion offers 
the most likely solution to eliminate dynamic overload- 
ing and to suppress the oscillations. 
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Letters to the Editor 


Dear Sir: 

I wish to commend Sears for bringing to light in his paper, 
“The Boundary Layer of Yawed Cylinders,” in the January, 
1948, issue of the JouRNAL OF THE AERONAUTICAL SCIENCES, the 
basic simplifications possible in the boundary-layer equations 
applicable to such bodies. In essence, the studies of Poritsky on 
yawed cylinders! have been completed by the examination of the 
boundary condition of zero slip and the establishment of the 
basic applicable Eqs. (2), (3), and (4). 

The heart of the paper would appear to be the recognition that 
“none of the physical quantities that describe the flow can be 
functions of y’’—i.e., span. Once that has been established, the 
findings that the boundary-layer profiles of the temperature and 
of the spanwise component of velocity are identical with that of 
the temperature during perpendicular flow conditions (surface 
temperature assumed constant and Prandtl Number assumed 
unity) automatically follow. Consequently, temperature bound- 
ary-layer profile solutions such as those of reference 2 describe 
also the side flow boundary-layer profiles for specific pressure 
gradients along yawed cylinder surfaces. Conversely, the side 
flow boundary-layer solutions given in the subject paper are also 
two-dimensional thermal boundary-layer solutions. The one 
shown for the region near the separation point—i.e., § = 0.7—is 
particularly interesting, showing as it does a temperature gradient 
and therefore heat transfer, at the separation point. Such a find- 
ing is in line with the results of references 2 and 3. On the other 
hand, we may borrow from the latter and state that the side flow 
boundary-layer profiles at different per cent chord stations will 
be similar when the pressure gradient along the yawed cylinder 
surface is favorable. 

Although the statement quoted at the beginning of the above 
paragraph seemed a likely assumption and is perhaps obvious,‘ 
I found it necessary, in order finally to convince myself of its 
validity, to construct a mental picture of a uniform ring airfoil 

of small thickness relative to the ring radius. Rotation of sucha 
ring around its axis while in a uniform flow field in the direction 
of that axis results in local flow conditions analogous to those on a 
yawed infinite cylinder, provided centrifugal force effects due to 
rotation are negligible. I find it self-evident that the local flow 
and boundary-layer conditions are identical all around the ring 
circumference. Since a yawed infinite cylinder corresponds to a 
rotating ring of infinite radius, it too must have uniform local 
flow conditions all along its equivalent ring circumference—that 
is, its span. 


These thoughts have been repeated not only possibly to help 
clarify the physical picture for others but also to point out that 
the boundary-layer results obtained on yawed cylinders are di- 
rectly applicable to spinning surfaces and to surfaces subjected 
to coaxial rotational flow fields when centrifugal effects are small. 
It would appear then that continued studies of the boundary 
layer on yawed cylinders will lead to clarification of the details of 
the boundary layer on such bodies as fixed cowlings and nacelles 
located behind tractor propellers, rotating cowlings, and rotating 
shells. 

It is a fortunate coincidence that Hitchcock’s paper, ‘‘Loss of 
Spin of Projectiles—Part II—Skin Friction Drag,’’ appears in 
the same issue of the JouRNAL; for in the light of the above re- 
marks, it is closely related to Sears’s paper. (Compressibility dis- 
tortions of the boundary layer: assumed effectively negligible in 
line with studies and experience to date in the Mach Number 
range dealt with.) Let us assemble the material to be used in 
studying Hitchcock’s paper: 


For THE LAMINAR BOUNDARY LAYER 


(a) The axial drag is independent of the spin speed—as indi- 
cated by Sears’s work. 

(b) In the absence of a pressure gradient, the local skin fric- 
tion drag, Dir, and the local spin-destroying couple, Mis, are re- 
lated through the equation 


Dir = 4Ui1Mis/di*w (I) 


where U; is the local axial velocity outside the boundary setts 
d, is the local diameter, and w is the spin speed. 

Eq. (I) is the analogy between the spin-destroying couple fie 
the skin friction drag corresponding to Reynolds Analogy be- 
tween heat transfer and skin friction when the Prandtl Number 
is unity and the surface temperature, constant. 

(c) When a pressure gradient exists, the local ‘‘momentum 
loss” drag, Diy, and the local spin-destroying couple, Mis, are 
related to a first approximation by Eq. (II). 


Dim = 4UiMis/di*w (II) 


Eq. (II) is the first approximation spin analogy corresponding 
to the momentum loss-heat transfer analogy in the presence of a 
pressure gradient.? 

In essence, (a) is the direct result of the uniform flow conditions 
in the spin direction, and (b) and (c) automatically result from 
the identity in the boundary-layer profiles of the temperature 
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(Prandtl Number, unity; wall temperature, constant) and the 
spin component of velocity. I should like to suggest that these 
conditions also exist in turbulent boundary layers. 

It may be readily shown that, in order for (a) and (b) to apply, 
the turbulent velocity and temperature fluctuations must be re- 
lated as in Eq. (III). 

Ui(Ts — Tw) U? Ui(wdi/2) 
Unfortunately there is little extant information clarifying the ap- 
plicability of (c). Je do know, however, that, from a theoretical 
standpoint, Reynolds Analogy, and therefore (b), is not correct 
when a pressure gradient exists.® 

Let us now examine the indications of experimental data. The 
following quote taken from page 162 of reference 6 confirms the 
applicability of (a): ‘“‘Ballisticians commonly assume the drag 
force is not affected by the spin. There are some measurements 
made on spinning models in German wind tunnels (during the 
war) which indicate that the effect of spin on drag is negligible 
for the spin rates produced by most service weapons.” 

The applicability of (b) cannot be so firmly established at pres- 
ent because of the lack of description of the detail geometry of 
the shells used in accumulating the data quoted by Hitchcock. It 
is perhaps significant, however, that the main conclusion of the 
paper—namely, that the overall skin friction drag coefficient de- 
termined by the loss of spin data depends on the shape of the pro- 
jectile, as well as the Reynolds Number—is predicted on the basis 
of the present considerations. For, neglecting pressure gradient 
effects, which fortunately occur most intensely only over the for- 
ward portion of the projectile where the surface area is small, the 
overall axial drag should be given by Eq. (IV) 


~ function of w (III) 


Dr = = (4/0) UiMis/di? (IV) 


Assuming a very slender projectile, the local axial velocity outside 

the boundary layer, U1, differs little from the free stream velocity, 
U, and Eq. (IV) becomes Eq. (IV)! 

Dr = (4 Mis/di (Iv)! 


The ratio of the summation on the right to the fraction, Ms/d?, 
used by Hitchcock in his Eq. (4), is purely a shape parameter. 
If Eq. (IV)! is correct, then Hitchcock’s breakdown of the data 
must show a scatter due to projectile shape wholly apart from 
pressure gradient effects. 


ADDITIONAL COMMENTS 


(1) It appears that the Reynolds Number used to define the 
skin friction drag should, strictly speaking, be based on the axial 
velocity, U, not the overall velocity, W, (see Eqs. (8) and (9) of 
Hitchcock’s paper). 

(2) The empirical factor, 0.45, used in Hitchcock’s Reynolds 
Number appears to account forthe average shape parameter dis- 
cussed above, as well as possibly a questionable definition of ef- 
fective fluid properties. Are the air density and viscosity values 
used those of the free air stream, those occurring at the surface 
of the shell, or those elsewhere? 

(3) The ballistician’s finding that spin does not affect the 
drag of projectiles tends to confirm the theory that the drag per- 
pendicular to the axis of a yawed cylinder having a turbulent 
boundary layer is a function only of the velocity component per- 
pendicular to the axis of the cylinder—exactly as in the case of the 
cylinder with a completely laminar boundary layer. 
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ARTHUR N. TirFoRD 
Section Engineer 
Aviation Gas Turbine Division 
Westinghouse Electric Corporation 


Dear Sir: 

The method of calculation of the suction force on the leading 
edge of a delta wing that has a “subsonic” type of leading edge 
which was presented in our paper, ‘‘Aerodynamic Performance of 
Delta Wings at Supersonic Speeds,” by A. E. Puckett and H. J. 
Stewart (JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 14, 
No. 10, p. 567, October, 1947) is, as was mentioned in the paper 
and by Leon Beskin (Letter to the Editor, JouRNAL OF THE 
AERONAUTICAL SCIENCES, Vol. 15, No. 3, p. 150, March, 1948), 
not entirely satisfactory. It is thus desirable that the conditions 
of applicability of this solution be discussed in a little more detail. 
Exactly the same limitations apply to the theory when it is 
applied to the leading edge of a sharp-nosed two-dimensional 
airfoil moving at subsonic speeds. As the correlation of experi- 
mental and theoretical results in this latter case is well known, it is 
best to discuss it first. 

As Mr. Beskin pointed out, the theory can hold at a mathe- 
matically sharp leading edge in a nonviscous fluid only if infinite 
suction pressures can be developed. On the other hand, a real 
wing must be expected to have a finite, but perhaps very small, 
radius of curvature. The linearized theory cannot, of course, 
be used to treat the details of the flow near the nose; however, it 
can be used to indicate the size of the region where nonlinear 
terms are important. If noncambered airfoils are considered, 
the flow at zero angle of attack is symmetrical, and there is a 
stagnation point at the leading edge. At a small angle of attack 
the stagnation point is shifted around on the lower surface of the 
leading edge, and the flow around the leading edge is, locally, a 
simple vortex flow. The low pressure over the surface produced 
by this vortex gives the suction predicted by the theory, provided 
the angle of attack is small enough so that the suction pressures 
are also much smaller than the static pressure. As the angle of 
attack increases, the suction increases, and the solution must 
break down by the time the theory (correcting for compressibility) 
shows zero absolute pressure on the surface. The maximum 
angle of attack for which the theory applies is thus determined 
by the leading edge radius; the greater the radius, the greater 
is the permissible angle of attack. 

From subsonic tests of sharp-nosed airfoils, it is found that the 
theory breaks down well before this zero pressure point because 
the large velocity gradients near the nose produce significant 
viscous forces. The viscous effects are generally observed to 
cause separation at the nose, and the separated layer then re- 
attaches itself to the wing surface a short distance back of the 
nose. The suction effect is practically lost in these circumstances. 

The same type of phenomena must be expected to occur for 
delta wings with large sweepback, O< k<1. The suction theory 
as presented should be valid for sufficiently small lift coefficients. 
The maximum angle of attack for which the theory is valid must 
be found experimentally. It is important to remember that 
rounding the leading edge should increase the maximum lift at 
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which this suction is available to reduce drag. Since it is possible 
that the optimum profile for minimizing zero lift form drag may 
also require a rounded leading edge, this idea may be of consider- 
able technical importance. 

H. J. SrEwart 


Associate Professor of Aeronautics and Meteorology 
California Institute of Technology 


Dear Sir: 

Congratulations are due George Gerard for his paper, ‘‘Opti- 
mum Number of Webs Required for a Multicell Box Under 
Bending” (JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 15, 
No. 1, p. 53, January, 1948). It is my own feeling that simplified 
approaches to optimum design have been insufficiently treated in 
the literature. 

It may be of interest to note how Gerard’s results would be af- 
fected by the use of a calculated value of k, to replace the as- 
sumed value of 4. N.A.C.A. Technical Note No. 1323, ‘Charts for 
the Minimum-Weight Design of Multiweb Wings in Bending,” 
which has become available since the completion of Gerard’s 
work, provides accurate values of k, for a variety of proportions. 
With these values, it is found that the optimum ¢,,/t, is equal to 
about 0.3),./b, for any values of b,./b, less than 2. Further, with 
the proportions so fixed, k, is accurately represented—again for 
bw/b, less than 2—by the empirical formula 


k, = + 0.03n*(b,./nb,)*] 
given here in terms of the number of bays m and the structural 


thickness ratio by/nbs. 
The solidity factor can now be written, using Gerard’s notation, 


which differs from the original expression by the addition of the 
first term in brackets and the substitution of the factor 0.3 for 
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0.4 within the second set of brackets. Preparation from this ex- 
pression of a chart (Fig. 1), like Gerard’s Fig. 2, shows that the 
range over which a given number of bays is the optimum is close 
to that previously shown. The results are given in the accom- 
panying figure. In no case does the number of bays prescribed 
for any structural thickness ratio differ by more than one accord- 
ing to the two methods. 

Designers will probably prefer the results obtained with &, as- 
sumed equal to 4, since this method prescribes either the same ora 
lesser number of bays than with k, taken as a variable. 


E. H. SCHUETTE 
Division of Tests 
Magnesium Laboratories 
The Dow Chemical Company 
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Inelastic Column Theories and an Analysis 
of Experimental Observations 


CHI-TEH WANG? 
New York Unwersity 


SUMMARY 


The inelastic column buckling theories as usually presented are 
found to be rather confusing. Many experimental results still 
cannot be explained satisfactorily by the theories as they stand. 
It appears to be pertinent, therefore, to define the inelastic buck- 
ling problem once more from a more rigorous mathematical point 
of view and to give a more rigorous mathematical treatment of 
the problem. The effect of assuming a constant tangent modulus 
on the buckling load in the tangent modulus formula is discussed, 
and many anomalies of column behavior in the inelastic region are 
explained. 

In making such a study, it is important to emphasize the dif- 
ference between an ideal column and an actual column and the 
difference between the buckling load and the ultimate load. 
Southwell’s method! for analyzing column tests, which was origi- 
nally proposed for the case of elastic buckling, is now shown to be 
valid for inelastic buckling. It is also shown that, in applying 
the method, instead of analyzing load and deflection measure- 
ments, simultaneous load and strain readings can be used. Thus, 
it is easier to measure the strain more accurately. 


INTRODUCTION 


M2 INTEREST IN INELASTIC COLUMN THEORY has 
recently been aroused by the presentation of 
Shanley’s paper.? The value of the paper lies in the 
explanation of certain physical facts concerning the 
inelastic buckling process. By means of an idealized 
hypothetical column, Shanley has shown that “‘bending 
begins at the tangent modulus load and the column 
load increases with increasing lateral deflection, ap- 
proaching the reduced-modulus load as a limit if the 
tangent modulus is assumed to remain constant.”’ 
There are many ambiguities in this statement; for 
example, what is meant by the phrase ‘bending begins” 
and what is meant by “‘the column load?’ And, still 
further, would such a conclusion, drawn from the study 
of an idealized hypothetical column, also be true for 
any other columns? 

A study of the experimental results on short col- 
umns—e.g., Fig. 8—immediately reveals the previous 
conclusion is not quite true, because, when the column 
compressive stress, P/A, is in the region of the pro- 
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portional limit of the material, the test results indicate 
that the so-called column loads are usually less than the 
tangent modulus load, while Shanley has shown that 
they must always be greater than, or equal to, the tan- 
gent modulus load. The experimental results also in- 
dicate that there is a consistent tendency for the column 
load to be close to the reduced or double modulus load 
rather than tangent modulus load when P/A is in the 
region of yield stress. There must be a satisfactory 
explanation for such consistent tendency, but it can- 
not be found in the present-day theories. 

It is the purpose of this paper to give a more rigorous 
study of the inelastic column buckling theories and the 
related experimental observations. By a more rigor- 
ous mathematical treatment, it is found that many 
anomalies of column behavior observed in the labora- 
tory can be satisfactorily explained. 


ELASTIC BUCKLING THEORY AND THE DIFFERENCE 
BETWEEN THE BUCKLING LOAD AND THE ULTIMATE 
LOAD OF AN IDEAL SLENDER COLUMN 


In order to clarify the discussions on inelastic column 
buckling, theories concerning elastic buckling will be 
briefly reviewed. 

The problem of elastic instability consists of deter- 
mining the smallest load at which the original straight 
form of equilibrium of a centrally loaded ideal column 
becomes unstable. To formulate the problem mathe- 
matically, it is usual to apply an infinitesimally small 
disturbance to the originally straight column and to in- 
vestigate whether this bent form of equilibrium can be 
maintained by the axial load, P, acting alone when the 
disturbance is removed. The idea of applying and 
removing the small disturbance, though of no impor- 
tance here, is important in the discussion of inelastic 
buckling load and will be elaborated to some extent 
later. Referring to Fig. 1, if such a bent form of 
equilibrium is possible, at any cross section perpen- 
dicular to the x-axis, the internal forces over the cross 
section can be reduced to a compressive force applied at 


the centroid of the cross section and a couple. The 
equilibrium equations are, therefore, as follows: 
s 
+ 
P 


Fic. 1. The bent form of equilibrium of an ideal column under 


the action of axial load P. 
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AE 
Fic. 2 (left). Theoretical load-deflection curves for ideal and 
actual slender columns. Fic. 3 (right). Bending of an ideal 
column under a lateral disturbance. 


JSacdA = P (1) 
and 
=0 or = —Py (2) 


where o; is the normal stress acting at an arbitrary point 
of the cross section which is at a distance y from the 
centroidal axis, A is the cross-sectional area, y is the 
deflection of the section from the centroidal axis, and 
the integration is taken over the whole cross-sectional 
area. Assuming that a plane perpendicular to the 
centroidal axis of the column remains plane after bend- 
ing, the x-component of the strain, ¢,, at an arbitrary 
point in a certain section is given by 


= — (n/R) (3) 


where ¢ is the strain at the centroidal axis, R is the 
radius of curvature, and 


1 (4) 
R [1 + (dy/dx)?)" 


where s is the distance along the centroidal axis and 6 
is the angle between s and x-axis. Integrating Eqs. 
(1) and (2) and denoting Ze by oo, since o, = Ee, (E 
being the Young’s modulus), it is found from Eq. (1) 
that.oo = P/A, and Eq. (2) becomes 


(d6/ds) + (P/EDy = 0 (5) 


where J is the moment of inertia of the cross-sectional 
area about the centroidal axis. The solution of Eq. 
(5) canbe expressed in terms of an elliptical integral 
and is given in any standard treatise on the Theory of 
Elasticity. When the deflection is small, the curvature 
can be approximated by d*y/dx?, and Eq. (5) is reduced 
to the well-known Euler’s equation as follows: 
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(d*y/dx*) + (P/ED)y = 0 (6) 


It is well known that the differential Eq. (6) will have 
a nontrivial solution only when 


P, = Cnr*EI/L?, n = 1, 2, 3, etc. (7) 
The smallest value of P is at nm = 1 and is equal to 
Po = CrEl/L? (8) 


where P,, is the so-called buckling load, L is the length 
of the column, and C is the coefficient of end fixity, 
Physically, it indicates that a bent form of equilibrium 
is possible when the axial load is of the magnitude given 
by Eq. (8). Solution of Eq. (5) reveals something 
more. It can be shown‘ that, at and below the buck- 
ling load, the column has only one form of equilibrium— 
i.e., the straight form. When the load is greater than 
the buckling load, there are two possible forms of 
equilibrium—the straight form which is unstable and 
the bent form which is stable. 

The results from Eqs. (5) and (6) are plotted in Fig. 
2. Suppose that the column is slender and that the 
material can sustain the compressive stress P/A = a 
without failure of elasticity. The loading of the ideal 
column will follow the curve OAB according to Eq. (6) 
and OAC according to the more exact Eq. (5). At some 
point D on the curve OAC, the center deflection 6 is so 
large that the sum of direct compressive stress and 
bending stress due to the bending moment Pé is suf- 
ficient to produce elastic breakdown. Beyond this 
point the actual curve OADE will begin to drop away 
from OAC, since the material can now sustain less load 
than before. The ultimate or maximum load a column 
can carry then occurs at point D, and the buckling load 
isat A. It is thus seen that for an ideal slender column 
the ultimate load is always greater than the buckling 
load. But aside from this fact, they are entirely dif- 
ferent theoretically, because the ultimate load is de- 
termined mainly by the mechanical properties of the 
material of the column, while the buckling load is de- 
termined mainly by the geometrical configuration of 
the column. The former is a load due to the material 
failure, and the latter is due to the instability of the 
straight equilibrium form. 


LOADING OF AN ACTUAL COLUMN 


Whereas an ideal column must be perfectly straight 
and made of material of homogeneous composition, an 
actual column is more or less imperfect in that it may 
be initially bent and may not be completely homogene- 
ous. It is obvious that the effect due to nonhomogene- 
ity of the material cannot be taken into account in the 
general theoretical discussion and will not be consid- 
ered in this paper. Assume that the column is not 
quite straight initially. Let + denote the initial de- 
flection of the column axis from the line of thrust. 
Then Eq. (6) is replaced by 


For ¢ 
discu 
is als 
will f 
ever, 
curva 
more 
ing of 
point 
colum 
curve 
Iti 
only 
is stal 
at wh 
the cu 
But t 
FG be 
mate 
Fig. 2, 
load f 
case O 
load is 
buckli 
propor 
able w 
the m 
part o 
AB. 
rence 
shorte 


ling 1q 


284 

both 

4; Prov 

1.€., 
\ 
4 
obta 
ic 
serie 
D 
1! 
| ly wher 
The 
; 
taine 
q 
Wh 


(6) 


INELASTIC COLUMN THEORIES 285 


[(d*y/dx*) — (d*yo/dx*)] + (P/E]|y = 0 (9) 


The form of y will now depend upon the form of yo, 
both quantities being regarded as functions of x. 
Provided that ‘yp vanishes at either end of the column— 
i.e., although the column has initial curvature it is cen- 
troidally loaded—a general solution of Eq. (9) may be 
obtained by expressing y and yo in terms of Fourier 
series. Thus, if 


y= 6, sin (nwx/L) (10) 


yo = 6, sin (nxx/L) 


where 6, and 3, are constants, we find on substitution 
that , 
5,,/[1 (P/P,)], (11) 


The deflection of the column at its center can be ob- 
tained by substituting x = L/2 in Eq. (10), or 


+a —... 


For an actual column corresponding to the ideal one 
discussed previously, the deflection versus load curve 
is also shown in Fig. 2. The loading of the column 
will follow the curve FG according to Eq. (12). How- 
ever, when the deflection is large, the simplified form of 
curvature is no longer a good approximation. If the 
more exact formula for the curvature is used, the load- 
ing of the column will follow the curve FIH. At some 
point J or J’ on the curve, the maximum stress in the 
column reaches the elastic limit, and then the actual 
curve FIJ.of FI’J’ will drop away from FH. 

It is to be noted that, for an actual column, there is 
only one form of equilibrium—i.e., the bent form which 
is stable. Therefore, there is no such equilibrium load 
at which an exchange of stabilities occurs. However, 
the curve FG does approach the line AB as an asymptote. 
But the loading curve will break away from the curve 
FG before the deflection becomes too large. The ulti- 
mate load, which corresponds to the point J or J’ in 
Fig. 2, can be either greater or smaller than the buckling 
load for the corresponding ideal column. Only. in the 
case of long columns can one expect that the ultimate 
load is near to the buckling load. This is because the 
buckling stress, P,,/A, in this case is well below the 
proportional limit and the deflection must be consider- 
able when the elastic limit is reached. Because of this, 
the maximum for P will occur somewhere on the flat 
part of the curve thereby close to the horizontal line 
AB. This agreement is a somewhat fortuitous occur- 
rence and cannot be regarded as a general rule for 
shorter columns. 


n = 1, 2, 3, etc. 


(12) 


INELASTIC BUCKLING THEORIES 


When P/A = opis beyond the elastic limit, the buck- 
ling load can be defined as the smallest axial load at 
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Fic. 4. A typical compressive stress-strain curve. 


which the bent form of equilibrium of an originally 
straight and centrally loaded column, resulting from 
the action of a small disturbance, becomes stable. The 
definition, that it is the smallest load at which the 
straight form of equilibrium of a centrally loaded ideal 
column becomes unstable, is not quite adequate. This 
can be seen from the fact that the column, once it is 
bent by a disturbance, will not return to its original 
straight form because of permanent set, even though 
the disturbance is removed. In other words, the 
straight form may be called unstable in such cases. 
In the elastic region, these two definitions are the 
same. 

As the understanding of inelastic buckling theories 
depends greatly on the idea of the application of an 
infinitesimally small disturbance, a word about such 
disturbance appears to be pertinent. It is to be noted 
that a disturbance is a small force or moment that can 
be applied and removed at will. Referring to Fig. 3, 
the disturbance is represented in the form of a small 
lateral load P’. If P’ is applied before P reaches the 
buckling load and is then removed, the ideal column 
will resume its original straight form when no fiber is 
stressed beyond the proportional limit. When some or 
all of the fibers are stressed beyond the proportional 
limit, the column, due to permanent set, will assume 
the bent form as indicated by the dotted line in the fig- 
ure, which is different from the original bent form when 
P’ is acting on the column. This indicates that the 
original bent form is not a stable form of equilibrium. 
If P’ is applied to the column when P is equal to the 
buckling load and then P’ is removed, the bent form 
will remain unchanged—i.e., this bent form is stable. 

While there is no importance in specifying when the 
disturbance is applied and when it is removed in the 
elastic case, the magnitude of the inelastic buckling 
load depends intimately on such actions, for the process 
of loading and unloading is now irreversible. .The final 
stress distribution at a cross section of the column can 
be the result of many different types of loading, de- 
pending on when the small disturbance is applied and 
when it is removed. The two extreme cases are as 
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Fic. 5. Two extreme types of strain distribution across column 
cross section. 
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follows: First, the column may be compressed to the 
critical stress, say, the point A on the stress-strain dia- 
gram (Fig. 4), and then the disturbance is applied and 
immediately removed. As a result, the column sud- 
denly bends. Because of the bending, the stress on 
the compression side is in the process of loading, fol- 
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lowing the curve AB, and the stress on the tension 
side is in the process of unloading, following the curve 
AD. The stress distribution in this case is shown in 
Fig. 5a. On the other hand, the disturbance may be ap- 
plied during the compression process before the com- 
pressive stress reaches the proportional limit and is re- 
moved when the load is increased to a value such that 
the bent form becomes stable. In this case, compres- 
sion and bending proceed simultaneously during the 
loading process, and there is no reversal of stress. The 
final stress distribution is shown in Fig. 5b. 

Denote E, = do/de the tangent modulus and o = 
P/A the stress at the neutral axis 0, Fig. 6a. Accord- 
ing to the first process of loading, von Karman’s double 
modulus formula can be derived. The derivation is 
well known but is briefly outlined here.for the purpose 
of discussion. 

Assuming that plane cross sections of the column 
remain plane during bending, the small bending stresses 
will be distributed along the depth of the cross section 
as shown in Fig. 6a. Since the bending is only slight 
in determining the buckling load, the stress curve on 
the loading side may be approximated by a straight 
line with a slope equal to E;. Eq. (1) becomes 


E So" ndA + E, Sie = 0 (13) 
and the bending moment is 


M = Jy omdA = (E/R) fo" n°dA + 
(E,/R) = (1/R)(El, + Ez) (14) 


in which J; and J; are the moments of inertia with re- 
spect to the neutral axis of the two portions of the cross 
section. (ZI; + E,J2) is usually written as E,J, where 
E, = + is the so-called von Karman’s 
double modulus. However, it is obvious that the ef- 
fect can be better described by writing (EI, + E,J2) = 
EI instead of E,J, since actually both E and J have 
been modified instead of changes in Z alone. Since 
the bending is only slight, EJ is approximately a con- 
stant throughout the length of the column. Eq. (2) 
then becomes 


(d*y/dx*) + (P/ET)y = 0 (15) 
and the buckling load is 


P., = Cr°EI/L? (16) 


This is essentially von K4rm4n’s double modulus for- 
mula except that E/ is now written in the place of E,J. 
In the case of colums with rectangular sections, it can 
be proved that (Fig. 6a), 


EI = [(4EE/(VE + VE) (17) 


Now let us consider the second process of loading. 
Let us again assume that the stress distribution over 
the section can be approximated by a straight line, in 
other words, we assume the tangent modulus is approxi- 
mately constant over the stress range and it can be 
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easily proved that the buckling load is given by the so- 
called tangent modulus formula as follows: 


= CrEI/L? 


where Z, is the tangent modulus corresponding to the 
stress P,,/A. The assumption of constant tangent 
modulus is usually justified by arguing that the bending 
is only slight; however, it would be interesting to see 
the influence of such an approximation on the magni- 
tude of the buckling load. A better approximation 
can be obtained by approximating the stress distribu- 
tion curve over the cross section by two straight lines 
instead of the usual one straight line, as shown in Fig. 
6b. Denoting £, and E, the tangent modulus corre- 
sponding to the stresses at two sides of the neutral 


(18) 


axis, the same results can be obtained as in the case of 
double modulus theory, merely by replacing E and E; 
by £; and EF», respectively. In the general case, the 
buckling load is equal to 


P., = CrEI/L* 


in which EJ = + 
tangular sections, 


EI = + 


(19) 


For columns with rec- 


(20) 


Write £, = ak, FE, = bE, where £, as before is the 
tangent modulus at o = P,,/A. For most engineering 
materials, a 2 1 and} < 1, and Eq. (20) can be written 
as 


EI = [4ab/(Va + = KEJ 
The factor 


(21) 


4(a/b) 
/at vob)? (a/b) +14+2Va/b 


is plotted against a/b for different values of } in Fig. 7. 
It is seen that, for low values of b which lie in the region 
of the “knee” of the stress-strain curve or near the 
proportional limit, the factor is likely to be less than 
one—i.e., the tangent modulus load is greater than the 
actual buckling load. For higher values of b which are 
in the region of large stress and strain or near the yield 
point, the factor is likely to be greater than one. The 
tangent modulus load is, therefore, less than the actual 
buckling load. This explains why the tangent modulus 
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theory predicts unconservative values when the buck- 
ling stress is near the knee of the stress-strain curves. It 
also is orfe reason why the tangent modulus theory pre- 
dicts values that are too low when the buckling stress is 
high (Fig. 8). 

Theoretically, for a centrally loaded ideal column, if a 
small disturbance is applied and immediately re- 
moved, the resulting bent form becomes stable when P 
is equal to the double modulus load, Py. Referring 
to Fig. 3, in this case, P’ is assumed to be applied and 
removed when P is equal to Pz. Whereas, if during the 
loading P’ is applied before P/A reaches the propor- 
tional limit and is removed when P is equal to P,,— 
i.e., load at which the bent form becomes stable—the 
buckling load, P.,, is close to the tangent modulus load, 
P,. If the disturbance, P’, is applied after P/A is 
beyond the proportional limit but below P,/A, then 
the buckling load would be equal to some value be- 
tween P, and P,. 

Since the tangent modulus, on which the values of 
EI depend, decreases rather rapidly shortly after the 
proportional limit is exceeded for most engineering ma- 
terials, the effect of plastic deformation on the load is 
larger than the effect due to the use of approximate ex- 
pression for the curvatures. The load-deflection curve 
would, therefore, drop immediately after buckling 
starts. It is thus seen that, for an ideal short column, 
the buckling load is equal to, or very close to, the ulti- 
mate load. 

If the testing of an ideal column is possible, there will 
always be disturbances in the testing machine during 
the loading process, and the double modulus load can 
never be obtained. As the magnitude of the disturb- 
ance created in a given testing machine is more or less 
constant, its effect on bending of columns depends ap- 
proximately upon the slenderness ratio L/p. For 
short columns with small L/p ratio, the column is rela- 
tively insensitive to small disturbances. In such 
cases, the column may bend when P/A is beyond the 
proportional limit and there is actually a reversal of 
stresses. The buckling load is, therefore, close to the 
double modulus load. For relatively slender columns, 
they are more sensitive to a disturbance of the same 
magnitude and will bend before the stress reaches the 
elastic limit, consequently with no stress reversal. In 
these cases, the buckling load is closer to the tangent 
modulus load. This effect can be seen clearly in Fig. 8. 


SOUTHWELL’S METHOD AND ITS EXTENSION 


Referring to Eq. (11), since P, = n*P,,, we have 


bn = 5,/[1 — (P/n*P..)) (22) 
As P approaches P,,, we see that 
& 4 9 
= = ‘i 


Therefore, 5; > 6: > 63 > ... 


By substituting in Eq. (12), it is evident that, if P is 
a fairly large fraction of P.,, the center deflection 6’ 
is approximately equal to 6,;—.e., 


= — (P/Pa)] (23) 


and Eq. (23) is a close approximation of Eq. (12), 
The P vs. 6’ curve approximates a rectangular hyper- 
bola having the horizontal line P = P,, as an asymp- 
tote. 

Since the deflections measured in testing are usually 
referred to the initial position, they are (6’ — 5) rather 
than 6’. Writing 6 = 8’ — 6, Eq. (23) can be rewritten 
as 


— & = 4,/[(P../P) — 1] (24) 
or 
P.(5/P) — 6 = (25) 


It is seen that if 5/P is plotted against 5, when P is near 
to P.,, the testing results will be approximately a 
straight line, and the inverse of the slope of this line is a 
measure of the buckling load P,,. Timoshenko’ also 
shows that this relationship holds true when there is 
some eccentricity in applying the load. This is the 
well-known Southwell’s method and has been widely 
used for estimating the elastic buckling load from test- 
ing results. However, because of the happy coinci- 
dence that the ultimate load of a column is close to its 
theoretical elastic buckling load, it is rather a common 
practice to regard the ultimate load as the buckling 
load, and the importance of Southwell’s method has 
not been properly emphasized. Actually, Southwell’s 
method represents a much greater achievement than is 
usually recognized. This is because, for imperfect 
columns, the buckling load is not defined, and all actual 
columns are more or less imperfect; thus, in the strict 
sense, actual testing of buckling is impossible. South- 
well’s method, however, provides a theoretically sound 
basis for analyzing the experimental data—from the 
test results of an imperfect column the buckling load of 
the corresponding perfect column can be estimated. 

In the case of inelastic buckling, Southwell indicated 
in his paper that his method would not apply, and it is 
so generally accepted. 

Let us refer to Eq. (9). If Pisa large portion of the 
buckling load P;, in the case of inelastic buckling, it can 
be replaced by the following relationship: 


[(d*y/dx*) — (d*yo/dx*)] + (P/E])y =0 (26) 


Assume that E/ is approximately a constant when P is 
close to P,,.. Then we have 


6 = — (P/P.,)] 
and 
— & = 4,/[(P../P) 1] (27) 
It is thus seen that the Southwell’s method can be ex- 
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Fic. 9. Calculated load-deflection curves for short columns 


(l/p = 75) with various amounts of initial deflection for mild 
steel with a yield point = 45,000 Ibs. per sq.in. (cf., von Karman‘). 


tended into the inelastic region, provided that the as- 
sumption of EJ being approximately constant can be 
justified. As EJ depends on the mechanical properties 
of the material, the justification of such an assumption 
should rely to a great extent on experimental evidence. 
The simplest way seems to be to apply the conclusion 
to the experimental results. If the experimental re- 
sults check with the theoretical prediction, it may be 


considered that the assumption is valid. This is done . 


by applying the extended Southwell’s method to the 
classical test results of von Karman,‘ the more recent 
ones of Gerard, and those of Horsfall and Sandorff ** 
It was found that the resultant plots are good straight 
lines, thus confirming the validity of the assump- 
tion. 

Donnell’? indicated that, even for columns buckled 
in the plastic range, the first part of the loading is al- 
ways elastic. By analyzing the measurements in the 
elastic range according to Southwell’s method, the 
buckling load can therefore be estimated. This argu- 
ment, however, is not satisfactory because the loading 
of a column in the elastic range is governed by Eq. (9), 
and the solution is given by Eqs. (10) and (11), where 
P., as obtained by Southwell’s method is Cr?EI/L? and 
not Cr*EI/L?. Also, if the buckling load is well in the 
plastic range, the readings in the elastic range may be 
those that are not the large portions of the buckling 
load. Although it will be shown later from the test 
results that the rectangular hyperbolic shape of P vs. 
6 curves is a good approximation, even well in the elastic 
range, this agreement can be regarded only as a some- 
what fortuitous occurrence. 

* When this paper was in preparation, the author received 


Horsfall and Sandorff’s report through the courtesy of F. R. 
Shanley. It was found that Southwell’s method was used in 


their report, but it was indicated there that the use ‘‘is a liberty 
not justified by theory.” 


It may be pointed out that, while the presence of ini- 
tial deflection in the elastic case does not materially in- 
fluence the ultimate load, the effect of initial deflection 
in the inelastic case is much more serious. This is be- 
cause, during the process of loading before the deflec- 
tion of the column becomes large, the effect of plastic 
deformation comes in to cause the ultimate load of an 
actual column to be lower than the theoretical buckling 
load of the corresponding ideal column. This was shown 
by von K4rman and is illustrated in Fig. 9. A method 
for estimating the buckling load from the experimental 
data in this case is even more important than in the 
elastic case, since the initial deflection of an actual 
column is usually difficult to determine, and, conse- 
quently, the accuracy of the test result cannot be 
properly estimated. The limitation of Southwell’s 
method in this case is thus evident. If the initial de- 
flection is too large, the experimental readings may be 
in the range where the loading curve is not approxi- 
mately a rectangular hyperbola. 


As Southwell originally proposed it, his method re- 
quires that the initial deflection reading be taken at 
zero load. In the vicinity of zero load, deflection read- 
ings are usually somewhat questionable. A more 
general method is suggested by Lundquist,* where the 
initial readings may be taken at any load less than the 
critical load. Lundquist proved that (6 — 4) vs. 
(6 — &)/(P — P)) is also a straight line, where (5 — 6:) 
is the amount by which the deflections are increased 
when the axial load on the column is increased from P 
to P. 1- 

Instead of measuring the center deflections, it is 
easier to measure accurately the strains by electric 
strain gages. If it is assumed that sections remain 
plane after bending, the curvature 1/R of a column at a 
given cross section is related to the difference in strain 
at two points, ¢ and &, on the particular cross section 
according to the equation 


Fic. 10. Modified ‘‘Southwell plot”’ of the test results of Schuette 
and Roy for a slender column. 
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TABLE 1 
A Comparison Between the Ultimate Load, Euler’s Buckling Load, and the Buckling Load Estimated by Southwell’s Method for 


Slender Columns 


Ratio of 
Slenderness Buckling Load Estimated Value 
Column Ratio. L Calculated Estimated Ultimate to Theoretical 
No. * Buckling Load* from Test Load, Kg. Value 
1 175.8 3,790 3,710 3,770 0.980 
2 146.0 5,475 5,453 5,430 0.995 
3a 116.2 8,645 8,590 8,630 0.994 
3b 116.1 8,610 8,758 8,750 1.017 
4a 103.0 10,980 11,220 11,160 1.022 
4b 103.5 10,920 11,090 0,860 1.015 
5 95.3 12,780 12,815 12,520 1.003 
6 91.3 13,980 13,750 13,580 0.984 


* Calculated from Euler’s formula by using the actual value of Young’s modulus measured by von Karman, i.e., 2,170,000 kg. per 


sq.cm. 


TABLE 2 


A Comparison Between the Ultimate Stress, Theoretical Buckling Stresses, and the Buckling Stress Estimated by Southwell’s Method 
for Short Columns 


Effective -——Calculated Buckling Stress-——— Buckling Stress 
Slenderness Tangent modulus Double moduius Estimated Uitimate 
Column No. Ratio, (L/p). theory theory from Test Stress 
von Karmén* Kg. per sq.cm. Kg. per sq.cm. Kg. per sq.cm. Kg. per sq.cm. 
Ta 88.1 2,400 2,690 2,780 2,760 
7b 88.0 2,400 2,690 2,780 2,685 
8 82.0 2,600 2,900 2,740 2,740 
9a 73.1 2,960 ,050 3,050 3,030 
9b 73.1 2,960 3,050 3,105 2,866 
10a 58.6 3,100 3,150 3,240 3,185 
10b - 58.6 3,100 3,150 3,130 3,080 
11 53.6 3,120 3,175 3,270 3,165 
12a 48.2 3,130 3,210 3,110 3,080 
12b 48.2 3,130 3,210 3,050 2,960 
13 47.3 3,140 3,215 3,100 3,060 
14b 38.2 3,160 3,320 3,480 3,320 
15a 28.8 3,220 3,560 3,700 3,395 
16 24.8 3,290 4,100 3,900 3 
17 22.0 3,450 Approx. 4,500 4,330 
,600 
Gerardf Lbs. per sq.in Lbs. per sq.in. Lbs. per sq.in. Lbs. per sq.in. 
1 21.7 41,900 47,500 44,900 44,880 
2 21.2 43,000 49,000 48,600 48,150 
Lbs. per sq.in Lbs. per sq.in. Lbs. per sq.in. Lbs. per sq.in. 
Horsfall and 
Sandorff f 29.9 37,000 42,000 37,600 37,200 


* Mild steel columns. 
24S-T aluminum-alloy columns. 


1/R = (4 — @)/t (@’/dx*)(y — yo) (28) 


' where ¢ is the width of the column across which the 


strain gages are attached. Differentiating Eq. (10) 
twice and combining it with Eqs. (24) and (28), the 
difference in strain measurements at the center of the 
column is 


Ae = @ — = — 1] (29) 


when P is near to P;. It is thus seen that, if Ae/P is 
plotted against Ae, the resultant will also be a straight 
line. This is confirmed by the test resultst of Schuette 
and Roy” in the elastic case (Fig. 10) and by the test 
results of Gerard in the inelastic case (Fig. 12). 

t These results were not given in the original report. The 


author is indebted to Dr. E, E, Lundquist of N.A.C.A. for his 
kindness in supplying the test data. 


APPLICATION OF EXPERIMENTAL RESULTS 


In order to test the proposed method of analysis, it 
is necessary to have related test values of load and cen- 
tral deflection for columns that have been loaded as 
centrally as possible. As was done by Southwell to 
test his method in the elastic range, von K4arm4n’s 
test results in the inelastic range are first to be analyzed. 

Von K4rm4n took special precautions to ensure ex- 
act centering of his applied loads, and for each column he 
tabulated his observations‘of load and deflections dur- 
ing the progress of the test. He classified his columns 
in three groups—described, respectively, as slender, 
medium, and thick. Slender columns are those having 
an L/p ratio greater than 90, and the buckling stresses 
of those are in the elastic range. Medium columns 
are those having an L/p between 45 and 90, correspond- 
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ing to those having their buckling stress between the 

portional limit and yield point. Thick columns are 
those for which L/p is less than 45 and the buckling 
stresses of which are greater than the yield stress of the 
material. p is the radius of gyration of the cross-sec- 
tional area of the column. 

The test data of the “‘slender’’ columns have been 
analyzed by Southwell. The buckling loads for the 
{deal columns thus estimated are listed in Table 1 as 
compared with the Euler’s theoretical buckling loads 
and the ultimate loads for the actual columns. 

The “medium” and “thick’’ groups of the test data 
are now to be analyzed. Instead of working with the 
buckling load, it was felt that working with the buck- 
ling stress was more convenient in these cases. Eq. 
(24) can be replaced by the following expression: 


5 = — 1] (30) 


where o = P/A ando,, = P,,/A. It is assumed that 
the cross-sectional area of the column remains constant 
with an increase of load. 6/o is plotted against 6 in 
Figs. lla and 11b. Except for columns 14a and 15b, 
where the recorded data are too few for such an analy- 
sis, it is seen that the test data lie closely on the cor- 
responding straight lines, thus confirming the applica- 
bility of the method. The buckling stresses so deter- 
mined are listed in Table 2. The buckling stresses cal- 
culated by the tangent modulus formula and the 
double modulus formula, as well as the ultimate stresses 
of the actual columns, are also tabulated in Table 2 
for comparison. The buckling stresses as computed 
by the double modulus formula were computed by von 
K4rman. The computation of the buckling stresses 
by the tangent formula is based on the average tangent 
modulus given by von Karman. It is seen that some 
of the buckling stresses so estimated are greater than 
the corresponding double modulus stresses. This is 
probably because the moduli used in the computation 
are the average values rather than the actual values for 
the particular specimens. A typical stress-strain curve 
tested by Gerard is shown in Fig. 13, and the variations 
of the moduli of a particular specimen to the average 
values can be large. 

In order to test the proposed procedure still further, 
Gerard, formerly of the Republic Aviation Corporation, 
has kindly supplied the author with the test data of two 
24S-T aluminum-alloy columns of rectangular cross 
sections (approximately 1'/, by '/2 in.). The columns 
were designed to fail in the plastic range and were manu- 
factured with great care to ensure as little initial curva- 
ture as possible. The specimens were tested flat- 
ended and were equipped with three electric resistance- 
type wire strain gages opposite each other on the wider 
sides, one pair being at the mid-section and the other 
two near the estimated inflection points. Referring 


to Eqs. (12) and (23), it is seen that, as load is in- 
creased, the column will bend approximately into a sine 
curve, because the first harmonic becomes large and 
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Fic. lla. ‘‘Southwell plot’ of von Karman’s test results for short 
columns. 
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Fic. 11b. ‘Southwell plot’? of von test results for 
short columns. 


004 of 
ae 
Fic. 12. ‘‘Southwell plot” of Gerard’s test results for short 
columns. 


other harmonics are little magnified by the load. Points 
of zero curvature can be approximately determined by 
passing a sine curve through Ae readings opposite each 
other along the length of the column. The effective 
length of the column is then the distance between these 
points of zero curvature or the inflection points. The 
test data are plotted in Fig. 12. It is seen that the test 
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Fic. 13. Compressive stress-strain curve of 24S-T aluminum 
alloy as tested by Gerard. 


40 
§ 
~ 06 i> 

04 

oz 
20 30 -60 .70 .60 
"70% in infle 


Fic. 14. Lundquist’s modified form of ‘‘Southwell plot’’ of the 
test results of Horsfall and Sandorff for a short column.® 


points lie closely on the respective straight lines. The 
buckling stresses so determined are listed in Table 2, 
together with the theoretical values and the ultimate 
stresses. 

Another set of accurately measured data is that of 
Horsfall and Sandorff,® of Lockheed Aircraft Corpora- 
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tion. Their data are plotted according to Southwell’s 
method modified by Lundquist. The result is repro- 
duced here as Fig. 14 and offers further evidence to the 
validity of the method. The test procedures are de- 
scribed in detail in references 2 and 5. 

It is interesting to note from the results in Table 2 
that all the buckling loads estimated by Southwell’s 
method are either equal to, or greater than, the corre- 
sponding ultimate loads, as they should be in these 
cases. 
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Static Thrust Analysis 


for Helicopter Rotors 


and Airplane Propellers 


WALTER CASTLES, JR.,* anp ARNOLD L. DUCOFFEt 
Georgia School of Technology 


SUMMARY 


This paper presents a solution for the relation between the 
blade angle, thrust, and torque coefficients of a helicopter rotor 
in hovering flight. This analysis presents the case wherein the 
induced velocity at a given radius, 7, is set up as a function of the 
thrust at that radius and a semiempirical constant, k, is intro- 
duced to account for the viscous shearing forces in the flow. 
Equations have been derived which yield, upon graphical inte- 
gration, the thrust and torque for the general case of blades with 
any solidity, twist, or taper. The solutions for the integrals of 
the equations for the special case of a constant chord, untwisted 
blade are also presented. The results apply to a helicopter oper- 
ating in or out of ground effect, and values of the semiempirical 
constant, k;, which have been calculated from experimental data, 
are presented for the entire range. This method gives values of 
the thrust and torque coefficients which are in good agreement 
with the available experimental data. Since the static thrust 
case of an airplane propeller is identical to that of a helicopter 
rotor in hovering flight, this method can be used for the calcula- 
tion of the static thrust and torque of modern variable pitch 
propellers. 


INTRODUCTION 


bee USE OF Glauert’s correction factor from his 
curve of 1/f vs. 1/F" is cumbersome, and the 
numerical ratios in hovering flight are of questionable 
accuracy because of the small-scale and wind-tunnel 
boundary corrections of the experiments from which 
they were determined. 

In an attempt to simplify the use of, and determine 
the position of, an experimental correction factor in the 
theoretical formulas, the following approach was 
evolved. 


NOTATION 


Angles in Radians: 

a, = angle of attack of a blade element measured between the 
zero lift chord line and relative wind 

x = angle between the inflow velocity at radius, r, and the 
normal to the rotor 

¢- = induced angle of attack of a blade element measured 
between the relative wind and the plane of rota- 
tion 


6 = blade angle for untwisted blades measured between the 
zero lift chord line and plane of rotation 
6, = blade angle at radius, r, measured between zero lift 


chord line and plane of rotation 
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Areas: 
dA = area of an elemental annulus, 2rrdr 
A = rotor disc area, rR? 


Forces: 
dL = lift on an elemental annulus 
dT = thrust on an elemental annulus 
T = rotor thrust 


Lengths: 
c = 
c= 
he= 
r= 


2 = 


Miscellaneous: 

slope of the lift curve of a blade element, per rad. 
mass flow through an elemental annulus 

angular velocity of the rotor, rad. per sec. 

mass density of air 


chord for a constant chord blade 

chord of a blade element at radius, r 

distance from rotor plane to the ground plane 
distance of a blade element from the axis of rotation 
distance from axis of rotation to blade tip 


a = 

aM’ = 

Q = 

p = 

Moments: 

dQ = torque about rotor axis of rotation of air forces on an 
elemental annulus 

Q = rotor torque about the axis of rotation 

Qs = profile rotor torque caused by the profile drag at a = 0 
about the axis of rotation 

Qs, = profile rotor torque about the axis of rotation caused by 
profile drag variation dependent on a 

Qs, = profile rotor torque about the axis of rotation caused 
by profile drag variation dependent on a? 


Q; = induced rotor torque about the axis of rotation 
Nondimensional Coefficients: 
Ca = rotor blade element drag coefficient 
Ca = profile drag coefficient of a blade element where 
2 
Con = 5, + (*) 
a a 
Ca; = induced drag coefficient of a blade element 
Ci, = lift coefficient on a blade element at radius, r 
Ce = rotor torque coefficient, Cg = Q/pxf?R5 
CQs, = rotor profile torque coefficient at a = 0 
CQs, = rotor profile torque coefficient dependent on a 
Cos, = rotor profile torque coefficient dependent on a? 
Cg; = rotor induced torque coefficient 
Cr = rotor thrust coefficient, Cr = T/prf?R* 
Nondimensional Parameters: 
b = number of blades 
ku *, = viscosity correction factor to account for the total 
° head change above the annulus (in and out of 
ground effect, respectively) 
ki = viscosity correction factor to account for the total 
head change below the annulus (out of ground 
effect) 
ku, ky = viscosity correction factors to account for the total 


head change below the annulus (in ground ef- 
fect) 


2 
hese 
erva- 
.oyal 
1932, 
\ero- 
irnal 
May, 
In- 
No. 
m of 
pp. 
Teh, ; 
947, 
Mc- 
No. 
i for 
pp. 
ctive 
time 
port 


ee 


AXIS OF ROTATION 


PV. 


02 ‘co 


Fic. 1. Hovering flow pattern—no ground effect. 


kis = ky + kn 
ky = viscosity correction factor in the final equations 
= rotor solidity for a constant chord blade, o = 
bc/xR 
Or = rotor solidity at radius, r, 0, = bC,/xr 
X = 0/Z 
Z = ak,o/16 


Pressures: 

total head 

atmospheric pressure 

static presure just above the rotor annulus 

static pressure just below the rotor annulus 

static pressure below the rotor annulus (in ground ef- 
fect) where the streamlines are normal to the rotor 


Velocities: 
V; = induced velocity at a blade element 
= component of V; normal to the rotor annulus 
V. = resultant velocity at a large distance below the rotor 
annulus where the streamlines are normal to the rotor 
plane 


V; = velocity below the rotor annulus (in ground effect) 
where the streamlines are normal to the rotor 
plane 

AV = total change in velocity along the line of thrust of the 


air in any elemental stream tube 


THEORY 


The mass flow through an elemental annulus of 
radius, 7, and width, dr, is pdAV,, where the normal 
component of the induced velocity, V,, through the 
annulus is V,; cos x (Fig. 1). Thus the mass flow 
through the annulus is 


dM’ = 2xpV; cos xr dr (1) 
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The thrust acting on the annulus will be the rate of 
change of momentum along the line of thrust of the air 
flowing through the annulus. In order to use Bernoul- 
li's equation to find the total change in velocity in 
terms of the induced velocity, V;, at the annulus, it is 
necessary to take into account the effects of viscosity, 
For the case of the helicopter in hovering flight at a 
very great distance from the ground, the velocity dis- 
tributions just below the rotor and at a great distance 
from the rotor are somewhat as shown in Fig. 2. The 
radial induced velocity distribution at the rotor is far 
from uniform, and well below the rotor it tends to be- 
come uniform. Bernoulli's equation applied to the 
flow in any stream tube assumes no loss in energy be- 
tween points (1) and (2) in Fig. 2, but in a real fluid 
there is a change in energy due to viscosity. As the 
velocity distribution is not uniform just below the 
rotor, there will be an exchange in energy between ad- 
jacent stream tubes and turbulence will result. Thus, 
Bernoulli’s equation, H = P + (*/2)pV?, will be modified 
by introducing a factor, k, to account for this change 
in energy. Then Bernoulli's equation can be written 
for any given stream tube as 


H = P + ('/2)kpV? (2) 
Writing the modified Bernoulli equation between free- 
stream conditions and a point just above the rotor 
(Fig. 1) gives 

P, = Py — ('/2)kypV? (3) 


Now writing the modified Bernoulli equation between a 


(1) 
ROTOR 


Vj 


(2) 


AXIS OF ROTATION 


Fic. 2. Radial induced velocity distribution at the rotor plane 
and a large distance below the rotor. 
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STATIC THRUST ANALYSIS FOR HELICOPTER ROTORS 


point just below the rotor and a point a very great dis- 
tance below the rotor (Fig. 1) where the static pressure 
has again reached free-stream conditions and the 
streamlines are parallel and normal to the rotor plane 
results in 


= Py + wo? — (4) 
Combining Eqs. (3) and (4) gives 
P, — Py = (*/2)kpV 2? (5) 
But 
dT = (P, — P,)dA (6) 
Also by momentum relations 
dT = dM’ AV = pdAV, cos (7) 
Combining Eqs. (6) and (7) gives 
V. = (2V, cos x)/k, (8) 
Thus 
dT = (2pdA V? cos? x)/k; (9) 


Making the usual approximations for small angles 
—i.e., sin g = tang = gandgos g = 1—the lift on the 
portions of the blades having chord, ¢,, lying within an 
annulus of radius, r, and width, dr, is, from blade ele- 
ment considerations, 


dL = (p/2)be,C;,2?r? dr (10) 


Eq. (10) contains the additional approximation that 
the rotational component of the resultant velocity at 
a blade element is negligible. Equating Eqs. (9) and 
(10) and letting o, = bc,/mr gives 


V; = (Qr/cos x) V (11) 


The induced angle of attack, ¢, at radius, 7, is, accord- 
ing to the previous approximation that tan g = ¢, 


g = (V, cos x)/Qr (12) 
Combining Eqs. (11) and (12) results in 
= Vo,C;,k,/8 (13) 


For the case of the rotor in ground effect, the stream- 
line pattern changes to one approximately as shown in 
Fig. 3. Writing the modified Bernoulli equation along 
a stream tube between a point in the free stream and a 
point just above the annulus gives 


Py = Po — (3/2) Rup 


Now writing the modified Bernoulli equation between a 
point just below the annulus and a point in the stream 
tube between the rotor and the ground plane where the 
streamlines are parallel and normal to the rotor plane 
results in 


Ps + = Ps + Ve" 


(14) 


(15) 


But P; can be written as 
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Fic. 3. Hovering flow pattern in ground effect. 


P3 = Po + (1/2) V3? (16) 


now letting k, = ky + Rx 


P, = Po + (1/2) Vs? (1/2) Rup Vi? (17) 
Combining Eqs. (14) and (17) gives 
P, — Py = (1/2)RypV3" (18) 


In a similar manner as shown in Eqs. (6) through (13), 
(for the case of the-rotor out of ground effect), the 
value of ¢ becomes 


= Va,C,k,/8 


The semiempirical factors, k, (out of ground effect), 
and k, (in ground effect), will assume different values 
for the stream tubes at each radius because of the non- 
uniform velocity distribution at the rotor. In this 
analysis k; and k, are replaced by the factor, , 
which represents a weighted average value of k, and 
k,, taken from r = 0 tor = R. Thus Eqs. (13) and 
(19) are replaced_by the following equation: 


(19) 


= (20) 
Since 
C, = aa, (21) 
then 
g = 0, — a, = 6, — (C;,,/a) (22) 
Combining Eqs. (20) and (22) gives 
6, — (C,/a) = (23) 
Solving Eq. (23) for C,,/a gives 
C, 320 ) 
1+ 24 
a 16 (24) 


The minus sign is used in front of the radical as the 
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STATIC THRUST ANALYSIS FOR HELICOPTER ROTORS 


angle of attack must be smaller than the blade angle of 
incidence. Thus, factoring Eq. (24) so that it may be 


plotted, 


16C, 

a*kyo, 
Figs. 4 and 5 represent a plot of 16C,,/a?k,o, vs. 160, + 
ako, from which values of 16C;,/a*k,o, can be obtained. 
By a previous approximation, the elemental thrust, 
dT, has been taken equal to the elemental lift, dL, thus 


T = So* dr (26) 
Eq. (26) reduces to 


Cr = 39 4 (7) 
Thus, with a given value of k, and given values of a, 
§,, and o, at each radial station, the thrust coefficient for 
the rotor with any twist or taper can be found by a 
graphical integration of Eq. (27). 

The rotor torque, making the same approximation on 
the induced angle, can be expressed as 


1 66, 
akyo, 


(25) 


Q = A® dr (28) 
The value of C, in Eq. (28) can be expressed as 
Ca = Ca + Cy (29) 
and 
Cu = (30) 


Combining Eqs. (28), (29), and (30) reduces to 
r\4 ( r 
3p,2 1 1 /s 4 

362 /o a*kyo, R R 
Values of C,, for the corresponding value of, C;, and 
Reynolds Number at each radial station can be ob- 
tained from the rotor blade airfoil profile drag polar for 
the airfoil at that station. 

Thus, with a given value of k; and values of Cy, a, 
s,, and 6,, at the various radial stations the value of 
Cg can be found by a graphical integration of Eq. (31). 

For the special case of the constant chord, untwisted 


blade, c, = ¢ = constant, o, = bc,/ar = o(R/r) and 
#, = @ = constant. Then, from Eq. (24) 


Letting 


Z=ake/16, X =0/Z 


then 


Also from Eq. (26) 
T = (0/2)bcC,,2*r? dr (34) 
Substituting Eq. (33) in Eq. (34) and integrating 
gives 
96Cr (1 — 3X)(1 + 2X)" 
X+5 5X? 


Thus, with values of 0, a, k:, and o for a given rotor with 
constant chord, untwisted blade, values of Cr can be 
found from Eq. (35). 

The rotor torque for the constant chord, untwisted 
blade is derived in two parts—induced torque and pro- 
file torque. 


(35) 


(a) Induced Torque 
The induced torque,-Q,, is defined as 
= So* dr (36) 


Replacing C,, by C, and substituting Eq. (33) for 
C, results in 


2,048Co, _ 

(1+ 2X) 2|- a- —3X)(14+2X) 1] (37) 


(b) Profile Torque 


The profile torque is defined as 
Qo = dr (38) 
Replacing C,, by: 
Cao = 50 + 51(Ci,/a) + (39) 
results in 
Coo = Com + Con + Cos (40) 


Combining Eqs. (38) and (39) and substituting Eq. (33) 
for C,, in Eq. (39) gives 


(1) Con = 0650/8 


32Ca, X 
(2) 6,ak\o? 4 + 3 


mike — 6X + 15X*)(1 + 2X)'* — 1| (42) 


(41) 


105X? 


2,048 Con yp , 16X 
++ + 4 


@ - 6X + 15X%)(1 + ox)" 2| + 


8 
105X? 


(1 — 3X) (1+ oxy” — 1] (43) 


el 


Thus, with values of a, 0, k;, and ¢ known, values of 
Cq can be found by the addition of Eqs. (37), (41), (42), 
and (43). 
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STATIC THRUST ANALYSIS FOR HELICOPTER ROTORS 


DISCUSSION AND RESULTS 


The numerical values for the semiempirical constant, 
h, for the rotor in ground effect* and out of ground ef- 
fect? were derived from experimental data. Fig. 6 
represents values of k, for a rotor at various distances 
from the ground plane. 

The calculation for the theoretical values of Cy; and 
Co were made assuming a value of the slope of the lift 
curve, @, of 27 per rad. Perfect fluid theory predicts a 
slope of 27 per rad. or greater® for a two-dimensional 
airfoil cascade, but the effects of viscosity tend to re- 
duce this value. Since the effect of viscosity has been 
included in the theory in the form of the factor, k,, the 
theoretical value of the slope of the lift curve of 27 per 
rad. is used. 

The blade airfoil profile drag coefficients for the rotors 
tested in the N.A.C.A. Full-Scale Tunnel* * were ob- 
tained from experimental data” * with the Reynolds 
Number based on the 3/,R. 

The tests conducted in the N.A.C.A. Full-Scale 
Tunnel* ° were run with rotors having radii of 12'/s ft., 
18 ft., and 19 ft. Assuming a slipstream contraction 
ratio of 2 to 1 it was found that the top of the balance 
house (30 ft. by 30 ft., approximately®), was larger than 
the minimum slipstream diameter. Accordingly, it was 
assumed that the ground plane was at the level of the 
top of the balance house, and values of k; were taken 
from Fig. 6 for the appropriate h/R ratio, where h for 
the Full-Scale Tunnel is 18.9 ft.® 

The theory derived previously was substantiated by 
comparison with experimental data?—* in the form of 
curves of Cy; vs. Cg (Figs. 7-11). The theoretical 
values are shown as solid lines, and the experimental 
data are plotted as points. The agreement between 
theory and experiment as shown in the figures appears 
satisfactory. 


CONCLUSIONS 


(1) The agreement between theory and experiment 
can be improved by the use of a semiempirical factor in 
the perfect fluid equations to correct for the effect of 
viscosity. 

(2) This theory predicts the thrust and torque for 
hovering flight and the propeller static thrust condition 
with a good degree of accuracy for any reasonable 
twist, taper, solidity, number of blades, and distance 
from the ground plane. 

(3) A value of the slope of the lift curve of 27 per 
rad. appears to be satisfactory for the case of hovering 
flight if the viscous effects are otherwise taken into ac- 
count. 


APPENDIX I 


Evaluation of 


As discussed previously the effect of viscosity is ac- 
counted for by the introduction of the factor, ki. As 
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shown in Fig. 6, the value of k, varies as the distance 
from the ground, h, is varied. In evaluating k; from 
the experimental data” * the slope, a@, was taken as 
2m per rad. For the case of the rotor in ground effect 
the value of h/R varied from 0.25 to 2.00? and for the 
rotor out of ground effect h/R was infinite.* 
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An Analytical Study of the Landing Shock 
Effect on an Elastic Airplane 


R. H. SCANLAN* 
Rensselaer Polytechnic Institute 


SUMMARY 


Several authors*: *: 7) 18 have considered the analytical problem 
proposed in this paper: given, the time history of the forcing 
functions (from theory’ or test experience®~" 12-14) applied to the 
main landing gear in vertical and drag directions during landing 
of an airplane. Assuming the airplane elastic, describe the vari- 
ous accelerations of its components. 

The solution presented here contains new features, some or 
all of which have been neglected by others. These are: 


(1) The inclusion of a sufficient number of degrees of free- 
dom adequately to describe important observed effects of actual 
landings. Such effects are exemplified by the “whipping” of 
landing gear, fuselage, and tail surfaces (occasionally of enough 
magnitude to cause faiiure). 

(2) The use of normal coupled vibration modes of airplane 
components as generalized coordinates. 

(3) The inclusion of the effects of structural damping through- 
out the complete time interval considered. 


The detailed motions of individual components in their proper 
phase relation to the other components of the structure are also 
considered. Of interest to the designer is the fact that over- 
conservative practices, such as the superposition of maximum 
effects, are avoided in the present analysis. 

The equations of motion developed here are fairly numerous. 
Since the forcing functions are not analytically prescribed, the 
solution of the equations of motion is lengthy, though straight- 
forward, by conventional approximate means. Use of a dif- 
ferential analyzer on the problem would, of course, be desirable, 
but the solutions are presented here in a way that is suited to ap- 
proximate calculation by hand methods. 


CooORDINATE SYSTEM AND KINEMATICS 


— A REFERENCE COORDINATE SYSTEM and other 
appropriate coordinates as given below. 


reference point at c.g. of airplane 

coordinate along centerline of airplane, posi- 
tive for x aft of R (in.) 

y = coordinate perpendicular to x-axis, axis parallel 
to ground during normal flight attitude, posi- 
tive outboard of R (in.) 

coordinate with axis vertical through R, posi- 
tive downward (in.). Fuselage vertical 
translation is given by 2 


R 
x 


The x, y axes are fixed in the airplane; the z-axis 
is fixed in space. . 
0; = fuselage rigid body rotation about R (rad.). 
Positive stalling 


Received November 7, 1947. 
* Associate Professor of Aeronautical Engineering. 


S;(x) = fuselage elastic bending deflection (in.), 
Positive down 
h,(x) = total vertical displacement of fuselage sec- 
tion at x (in.). Positive down 
h.(y) = total vertical displacement of wing section 
y (in.). Positive down 
6.(y) = total rotation of wing section y (rad.), 
Positive stalling 
Sw(y) = vertical bending component of wing elastic 
deflection along y-axis (in.). Positive 
down ; 
torsional component of wing elastic deflec- 
tion along y-axis (rad.). Positive stalling 
Yo = point along wing where main landing gear is 
located (in.) 
r = distance from wing downward along landing- 
gear strut (in.). Positive down 
tip = Value of r at wheel axle 
h,(r) = horizontal travel of section r of oleo strut 
(in.). Positive forward 
s,(r) = horizontal elastic deflection of section r of 
oleo strut (in.). Positive forward 
s;(y) = elastic vertical bending deflection of stabi- 
lizer “elastic axis’? at section y (in.). 
Positive down 
total vertical travel of stabilizer ‘‘elastic 
axis” at section y (in.). Positive down 
-. = x-coordinate of stabilizer ‘‘elastic axis” 
(in.) 
The airplane will be assumed initially stationary in 
space, with rigid body motions of translation (z) and 
rotation (6;) possible. The fuselage will be assumed 


elastic and free to bend vertically. The wing will be where 


assumed elastic and free to bend and twist. The 
landing gear will be assumed an elastic beam cantilever 
from the wing. The horizontal stabilizer will be as 
sumed an elastic beam free to bend but not to twist. 
Under the above assumptions the following relation 
hold: 


hy (x) = 2 + + 

holy) = 2 + Sy (y) 

Bw(y) = Oy + (1 
h,(r) = + 

hay) = 2 + + + 5,(y) 
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KINETIC ENERGY AND GENERALIZED COORDINATES 


The kinetic energy of the entire airplane system may 
be expressed as the total of various components. The 


result is: 


T='/s (x)dx + Sy 4 s x76, (x)dx 
+ s + Sp (x)dx + 
SG 5 my (x)dx + Sp + 
5,’m,(y)dy +S s,m,(y)dy + 4 (2) 

Sot (hw)? Mo(y)dy + So+ ghubwSo(y)dy + 

1/9 5,?m,(r)dr } 


where 


J; +s indicates integration throughout fuselage and 
stabilizer as though stabilizer were rigidly part of 
the fuselage 

J; indicates integration throughout stabilizer alone 

Jo+,g indicates integration throughout wing and 
landing gear as though landing gear were rigidly 
part of the wing 

J; indicates integration throughout landing gear 
alone 

m,(x) is fuselage-stabilizer mass at section x (Ibs. 
sec.?/in.?) 

m,(y) is stabilizer mass at section y (Ibs. sec.?/in.?) 

M,,(y) is mass of wing section y (Ibs. sec.?/in.?) 

Ss(y) is wing static unbalance about y-axis at section 
y (positive for c.g. aft of y-axis) (Ibs. sec.?/in.?) 

I, (y) is wing mass moment of inertia per unit span at 
section y (Ibs. sec.”) 

m,(r) is landing-gear mass per unit length (Ibs. sec.*/ 
in.”) 


Let generalized coordinates be introduced by the rela- 
tions 
Sw = Inti + Inks 
Oy = aki + 
= hsk3 
Se = (3) 
5 = 
2 = & 
6; = J 
where 


hy, he are nondimensional functions of y proportional 
to the bending displacement components at the 
y-axis of the first and second cantilevered wing 
coupled vibration modes, respectively. These 
wing modes are calculated as from a rigid canti- 
lever support at the root and are referred to point 
R. hh = ky = 1 at the wing tip. 

@, a are functions of y (rad. per in.) proportional to 
the torsional displacement component about the 
reference axis of the first and second cantilever 
wing coupled modes, respectively (see /;, hz, above). 
a = a, = 1 at the wing tip. 


hs is a dimensionless function of x giving the un- 
coupled fuselage first bending vibration mode, rela- 
tive to point R. For this mode the fuselage is as- 
sumed independent of wing and stabilizer, other 
appendages being assumed rigid and integral with 
the fuselage. hs; = 1 at the trailing edge of the 
fuselage. 

hy is a dimensionless function of y defining the dis- 
placement in the first uncoupled bending mode 
along the “elastic axis’ of the stabilizer. This 
mode is calculated assuming the stabilizer canti- 
lever from a rigid support. /, = 1 at the stabilizer 
tip. 

hg is the dimensionless function of r giving displace- 
ment in the first horizontal bending vibration 
mode of the oleo strut. hs = 1 at the wheel axle. 

&, & generalized coordinates, functions of time ¢, 
defining wing coupled mode motion. 

£3, £4, generalized coordinates, functions of time 
representing fuselage bending, stabilizer bending, 
and landing-gear bending, respectively. 

generalized coordinates, functions of time ¢, 
representing airplane free body translation and 
pitch, respectively. 


The result of introducing these generalized coordi- 
nates is that 


hy = ig + + hake 
6 


h 6 + xtr + Nsks (4) 
hy = + + + 
h, = + + As(xs)és + Maks 


The kinetic energy, Eq. (2), can now be written 
7 7 ? 
t=lj= 


where Ay = Ax and 


An = Jorg + 2 So(y)dy + 
earl o(y)dy 

An = dove + 2 bate + 

Ay = +5 hs?m (x)dx 

Au = JS; hem,(y)dy 


) 


As = hs?m,(r)dr 
Acs = M, = mass of entire airplane [(Ibs./in.) 
sec.?] 


An = In, = mass moment of inertia of entire } (6) 
airplane about reference axis (Ibs. in. 
sec.”) 

= 0 (due to orthogonality of normal modes) 

Ai = 0 (no coupling assumed between wing 

modes and fuselage bending) 

Ay = 0 (no coupling assumed between wing 
modes and stabilizer bending) 

Ars = an(y,) S; rhem,(r)dr 

= Sor, himey)dy + Sor, 
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= Sore + aale(y)dy ) 

Ass = 0 (no coupling assumed between wing 
modes and fuselage bending) 

Am = 0 (no coupling assumed between wing 
modes and stabilizer bending) 

As = rhym,(r)dr 

= + a2Sy(y)dy 

= + Sos 

Ay = hs(x,) ham.(y)dy 

Az, = 0 (no coupling assumed between fuselage 
vertical bending and landing-gear hori- 
zontal bending) 


Ax = (x)dx 


= xhym (x)dx > (6) 

As = 0 (no coupling assumed between stabilizer 
vertical bending and landing-gear hori- 
zontal bending) 

Ag = JS, hum,(y)dy 

Ag = hyn,(y)dy 

Ass = 0 (no coupling assumed between transla- 


tion and landing-gear horizontal bend- 
ing) 

As; = J; rhsm,(r)dr 

Ag; = 0 (no coupling assumed between free body 
translation and free body pitch—i.e., 
plane is in static equilibrium under 
sufficient external forces 


POTENTIAL ENERGY 


This is a quadratic function of those displacements 
that store elastic energy. If two normal modes are 
orthogonal, no cross product terms are present in their 
contribution. Rigid body motions contribute no po- 
tential energy. As a result, 


5 
U='/, Bé? (7) 
‘= 
where, with recourse to the free vibratory cases where 
§,(j = 7) are all constrained, the B, may be evaluated as 
B, = 
the value w, being the natural (undamped) frequency 
associated with the mode &;. 
GENERALIZED FORCES 


Consider, at a single landing-gear wheel, two exter- 
nally applied forces that may be considered known 
functions of time applied at the axle of the wheel. 
These are: 


D(t)—‘‘spinup”’ force horizontal (positive forward) (8) 
V(t)—oleo reaction vertical (positive down) 


It is sufficient to consider only one-half the entire air- 
plane throughout the analysis, merely selecting appro- 
priate vibration modes and forcing functions to ac- 
company symmetrical or unsymmetrical landing con- 
ditions. 


The generalized forces Q; are given by 


Q:; = Dit) (x {reiplér + + a2(y,)&] + 
) + Vit) (= {& + + + 


+ [aly,) + ) 0) 


and their specific values are 

Q1 = [reipor(y,)] + [n(y,) + 

= 0 


Q=0 (10) 
Qs = D(t) [hs(rrip)] 

Qs = Vit) 

Q; = [D®]nip + [VO] x, 


DAMPING AND THE EQUATIONS OF MOTION 


Several authors‘ 16 have discussed internal 
damping in elastic structures. It isnot viscous; there- 
fore, another type will be introduced here. Theodorsen 
and Garrick,'! following Becker and Féppl, employed a 
damping force proportional to displacement and in 
phase with velocity of a harmonically oscillating struc- 
ture. Thus, they introduced a damping force 


P, = +igmoth (11) 


where mw*h is inertial force. Reissner‘: * has suggested 
a means of introducing the same damping coefficient, g, 
without use of complex quantities. His suggestion 
amounts to replacing w* in the typical undamped har- 
monic equation 


x + wx = 0 
by 
w*[1 + (g/w)(0/dt)] 


This technique will be employed here. 
Application of the Lagrangian Equations in the form 


(d/dt)(OT/0é,) + (0U/0&) = 
to Eqs. (5), (7) and (10) yields the following equations: 
+ + Arsés + Arsés + = Qi(t) (12a) 
Asst: + rt, + Aasts + + Ant: = Q2(t) (12b) 


+ ws*Aasts + Assts + Assts + = 0 = Q,(t) 
(12c) 


Auks + wPAuks + Asks + Assts + Aut: = 0 = Qult) 
(124) 


Assis + ws*Assts + + + Asrtr = Qs(t) (124) 
Agsts + + Arste + Assts + Asst: = (12f) 


+ + And: + Aséts + Ants + Asrts = Qr(0) 
(12g) 
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LANDING SHOCK EFFECT ON AN ELASTIC AIRPLANE 


Here damping has not been introduced. It is found 
simpler to introduce it at a later stage in the solution of 
these equations. It may be pointed out that the 
later techniques used (in particular, introduction of 
normal coordinates) are not necessary for solution if 
Eqs. (12) are solved mechanically by means of a dif- 
ferential analyzer. If this is done, it is feasible to in- 
troduce damping terms of the form gw, directly into 
each of equationsz = 1, ...., 5 in Eqs. (12) before solu- 
tion. These terms are omitted here. Their introduc- 
tion would make solution of Eqs. (12) more compli- 
cated. Instead, the damping will be associated with 
the normal coordinates to be employed later in solu- 
tion. When this is done the reasonable assumption is 
implicitly made that the damping effect that serves to 
couple normal coordinates (natural modes) with each 
other is negligibly small in comparison with the other 
disturbances that excite the normal modes. 


SOLUTION OF DYNAMIC EQUATIONS 


Eqs. (12) may be easily reduced to five equations 
with %, = removed. Replacing Eqs. (12) by a single 
equivalent matrix equation in {, ...., yields 


AE + BE = X(t) (13) 


where the 5 by 5 matrices A and B and the single column 
matrix X(t) will have obvious definition. If K is the 
square matrix of the modal columns of C = BoA, 
then? 

KOCK =A 
where A is the diagonal matrix whose nonzero elements 
are the characteristic roots \; = 1/0,? of the system 


(13) and the first five characteristic roots of the system 
(12). Thus, introducing 


t=Kx (14) 
results in defining normal coordinates x of the system, 
Eq. (13) becoming: 

E+ AE = = (15) 
Solution of Eq. (15), without damping, as it stands, 
and later modified to include damping, will occupy our 


attention. 
Equation 7, from Eq. (15), may be written 


Zi + = Y,(é) (a 5) 


where matrix notation is now dropped. 
This may be solved,’ without damping, by 


= (1/a,) sin — 7)}Vi(r)dr (16) 
with the results, by differentiation: 
= cos [a(t — dr (17) 


= — & sin — 7) Vi(r)dr + Y(t) (18) 
The matrix Eq. (14) yields 
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(19) 


from which accelerations in the coordinates £ can be ob- 
tained. 

If it is now desired to apply Eq. (15) to a damped 
structure, let g; be the damping coefficient associated 
with the ith normal mode. (This may be assumed in 
practice to be close to the value g; associated with 
the generalized coordinate £,.) Then the typical 
equation of Eq. (15) is, as indicated earlier, represented 
by 


Ki 


[1 + = 
The solution of this* is given explicitly by 
t 
sin [a, W1 — (g2/4) (t — ¥(r)dr (20) 
and, as before, velocity and acceleration are given by: 


i, = £ [V1 = | 


0 


sin @V 1 — (g,?/4)(¢ — Yu(r)dr 


0 


[a V1 — (2/4) + 2V1 — (¢/4) 


sin 1 — (g,7/4) — + 
gw, COS V1 — (g,?/4) (t — Vilr)dr | 


The values of §, are again given by Eq. (19). For g, 
= i vi — (g:?7/4) = 1 may be used to simplify 
results. 

Comparison of damped with undamped results above 
shows clearly the frequency, phase, and amplitude-ef- 
fects of damping as assumed present. 
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Afterburners for Turbojet Engines’ 


F. A. CLEVELAND?{ 
Lockheed Aircraft Corporation 


ABSTRACT 


The general problem of pure jet aircraft take-off distance is in- 
dicated, and the basic reasons for the excess-thrust deficiency as 
compared to propeller aircraft are shown. The inherent advan- 
tages and disadvantages of several methods of short-duration 
thrust augmentation, such as fluid injection or rocket boosters, 
to alleviate the difficulties are mentioned. The afterburner con- 
cept is introduced and discussed; thrust augmentations of 45 
per cent low speed and over 100 per cent high speed are shown to 
be feasible. A method for analyzing afterburners is presented, 
including the thermodynamic derivations required. An example 
comparative performance analysis of a hypothetical jet-propelled 
aircraft, with and without afterburners, shows sea-level rate-of- 
climb increases of 200 per cent maximum speed increases of 10 
per cent, and normal take-off distance reduction of 50 per cent 
with afterburning. The example concludes the discussion. 


INTRODUCTION 


NE OF THE MOST SEVERE PROBLEMS encountered in 
the design of turbojet-powered aircraft is that of 
take-off within reasonable distances at fairly heavy 
weights. Unfortunately, the one characteristic that 
makes the pure jet power plant so well suited for high- 
speed aircraft—its nearly constant thrust at all air 
speeds at a given altitude and engine speed—lmits 
sharply the aircraft’s acceleration during ground run 
and transition climb. Though the thrust available 
from a propeller power plant falls off at very high air- 
plane speeds, because of compressibility effects, it 
doubles or triples its normal flight value at a given 
engine brake horsepower when the airplane velocity ap- 
proaches zero—that is, during the take-off procedure. 
In addition to the converse characteristics of the 
two types of power plants indicated, the turbojet is 
designed to give optimum fuel consumption at near 
maximum thrust conditions, and thus the airplane in 
which it is installed is designed for high-speed cruising. 
The take-off thrust available is generally not more 
than 15 per cent greater than the normal rated thrust. 
Reciprocating engines, on the other hand, cruise at 
only about 50 per cent take-off power, so that an air- 
plane designed for this cruise power has a large margin 
of excess power available for take-off. The turbojet 
suffers a further penalty in take-off operation in that 
the air inlet duct entrances are designed for high-speed 
conditions; in the low-speed take-off situation, stall- 
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ing of the duct lips causes marked decreases in the ram 
ratio available to the engine, with accompanying severe 
thrust declination. 


Some SOLUTIONS OF THE AUGMENTATION PROBLEM 


Several means of overcoming this inherent problem 
in turbojets have been proposed and tested with vary- 
ing degrees of success. In general, it may be stated 
that devices for take-off thrust augmentation also may 
be a great asset to military aircraft in that they provide 
increases in high speed and climb performance for 
emergency use. 

Auxiliary turbojet units may be installed in the air- 
plane. Such an arrangement provides the extra thrust 
desired but adds considerably to the equipped weight 
of the airplane and, when not in operation, imposes drag 
increments deleterious to high-speed aircraft. Powder 
rocket booster units serve admirably in- take-off, but 
they are not at present controllable after ignition, are 
expensive for continual use, and do not provide emer- 
gency high-speed augmentation. Liquid rockets will 
probably solve the controllability and emergency speed 
difficulties, but they introduce the weight and space 
problem and the difficulty of carrying several auxiliary 
fuels. 

Probably the most satisfactory solution flight-tested 
to date is the injection of a fluid spray into the turbojet 
compressor inlet or outlet. Extensive test data indi- 
cate that up to 25 per cent, increases of normal thrust 
are obtainable with moderate quantities of water in- 
jection. However, this method of augmentation is 
limited by combustion chamber flame blowout at a 
certain maximum liquid flow rate. It requires liquid 
tanks separate from the engine fuel supply and appar- 
ently decreases in usefulness as the airplane gains alti- 
tude. 

A fundamentally different approach to the problem 
is obtained with the notion of the turbojet afterburner. 
It is believed that this idea overcomes partially or 
wholly all of the disadvantages outlined heretofore for 
other methods of augmentation, and the progress cur- 
rently being made in several research and development 
laboratories indicates impending completion of practical 
units. 


THE AFTERBURNER CONCEPT 


Development of the Idea 
The basic limitation of the thrust output of a given 
turbojet engine is the maximum temperature permis- 
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Fic. 1. Afterburner configuration. 


sible at the turbine inlet. This defines the maximum 
stagnation temperature rise available from the compres- 
sor inlet to the nozzle exit, and it is this temperature 
rise, combined with the mass gas flow, which deter- 
mines roughly the maximum thrust of the engine. In 
order to keep the turbine inlet temperature below its 
limit (about 2,200°R. in present practice) in normal 
engine operation, 200 to 400 per cent excess air over that 
required for combustion is pumped through the unit. 
If a higher temperature increment is to be made some- 
how attainable, it is apparent that this air could be 
utilized as the oxidizer of additional fuel. 

The development following this idea is that of the 
afterburner, a combustion chamber aft of the turbine 
(or turbines) in which further burning and heating may 
take place. Effectively, this additional burner is a ram- 
jet, utilizing the high-pressure, high-temperature ex- 
haust gases from the turbojet. Since such a unit has 
no moving parts and thus imposes extremely small 
loads on components exposed to the hot gases, practi- 
cally the only temperature limit encountered is the ulti- 
mate red-heat temperature of the afterburner shell. 
And since this shell may be cooled to a certain extent 
by several means, the maximum allowable gas tempera- 
ture in the burner may be in the neighborhood of 4,000° 
R. Under such conditions, thrust increases between 80 
and 120 per cent are possible at airplane speeds on the 
order of 600 m.p.h., and increases of nearly 50 per cent 
accrue for the critical take-off condition. 

A schematic drawing of one type of afterburner is 
shown in Fig. 1. Its configuration is essentially that 
of a simple ram-jet: a diffuser section mounting fuel 
spray bars or nozzles; an open-grid flame arrestor, 
followed by a cylindrical combustion chamber; and a 
variable-area exit nozzle. The turbojet exhaust gases 
slow down in the diffuser, converting some of their 
kinetic energy into static pressure and simultaneously 
vaporizing and mixing with the fuel spray from the 
nozzle bars. Steady-state combustion is initiated 


downstream of the flame arrestor, resulting in a high- 
tea perature, high-velocity jet at the exhaust nozzle 
exit. 
possible in the unit. 
injection at several different stations; 


It should be noted that many configurations are 
Prevaporized fuel injection; 
direct nozzle 


Subscripts 
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burning without a flame arrestor; and myriad combj- 
nations of different diffusers, flame arrestors, and com- 
bustion chambers all have various unique merits that 
recommend them for specific installations. The con- 
siderations and statements made in the analysis to fol- 
low are thought to be of general enough nature to ap. 
ply to most of them. 


Inherent Advantages and Disadvantages 


The afterburner itself is relatively light; preliminary 
designs indicate a weight equal to about 20 per cent of 
the weight of the turbojet engine, which it supplements, 
Extremely small drag increases are encountered for 
most installations, since the afterburner is merely a co- 
axial extension of the jet engine, generally of equal or 
smaller diameter. The amount of augmentation ob- 
tained may be varied between wide limits at the will of 
the operator and is, so far as is known, available at any 
aircraft speed or altitude at which the turbojet will 
function. Since the afterburner and turbojet fuels are 
identical, no separate tanks are required for the former. 

The specific fuel consumption of the afterburner is 
considerably higher than that of the turbojet—values 
between 2.0 and 3.0 lbs. per hour per pound of thrust 
are indicated for the afterburner-turbojet combination 
—but its use is intended only for short periods when 
high performance is required. Some losses occur in the 
exhaust gas stream from the turbojet when the after- 
burner is not operating, but the thrust decrease is rela- 
tively small, less than 5 per cent. 


NOMENCLATURE 


duct cross-sectional area, sq. ft. 

local speed of sound, ft. per sec. 

specific heat for average static temperature, B.t.u. 
per lb. per °F. 

jet thrust (gross thrust), Ibs. 

net thrust, Ibs. 

acceleration of gravity, ft. per sec.” 

lower heating value of fuel, B.t.u. per 1b. 

Mach Number. 

mass flow, slug per sec. 

total pressure, lb. per sq.ft. 

static pressure, Ib. per sq.ft. 

gas constant, 53.3 ft.Ib. per Ib. °F. 

total temperature, °R. 

static temperature, °R. 

velocity, ft. per sec. 

air flow, lb. per sec. 

fuel flow, lb. per sec. 

gas flow, lb. per sec. 

gas state parameter 

gas velocity parameter 

specific heat ratio 

combustion efficiency 

mass density, slug per cu.ft. 


gas jet state when expanded to ambient static pres- 
sure 

station immediately aft of turbine 

station at combustion chamber entrance 

station at exhaust nozzle exit 
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AFTERBURNERS FOR 


0 = ambient static station 

TJ. = turbojet 

AB. = afterburner 

Tot. = combined turbojet and afterburner 


ANALYSIS OF AN AFTERBURNER 


Outline of the Approach 


This entire analysis is intended to apply to the situa- 
tion generally encountered in afterburner design—that 
is, when the basic performance data of a specific turbo- 
jet engine are known and the overall performance of a 
proposed afterburner design to be added is desired. 
Included are the several subsolutions for temperature, 
pressure, and velocity at various points which may dic- 
tate design alteration. This technique is, of course, 
equally applicable to any hypothetical turbojet for 
which complete performance has been assumed or calcu- 
lated. 

The turbojet performance data furnished by the en- 
gine manufacturer apply to a particular engine when it 
is equipped with a specified exhaust nozzle or tailpipe 
and nozzle. In order to simplify the subsequent no- 
menclature, the stations in the standard tailpipe and 
nozzle are given the same numerical designation as 
their counterparts in the afterburner, but primes are 
added. Thus the turbine-outlet-station total pressure 
is P; when the afterburner is in place and P;’ when the 
engine tailpipe is connected; similarly, with the exhaust 
nozzle outlet, station 10. Reference is made to Fig. 1 
for station numbers. 

The analysis is made in seven basic steps; they are 
listed briefly here to offer better perspective of the de- 
tails to follow: 

(1) The exhaust gas state conditions at the stand- 
ard tailpipe nozzle exit, station 10’, are calculated from 
the performance data on the engine at a given altitude, 
air speed, and engine r.p.m. 

(2) Total pressure losses are calculated or assumed 
between stations 5’ and 10’, and the gas state condi- 
tions at 5’ are determined. 

(3) The conditions calculated at 5’ are established 
as the conditions at 5 (i.e., when the standard tailpipe is 
removed and the afterburner substituted). 

(4) Total pressure losses are assumed or calculated 
between 5 and 6 and between 6 and 10. 

(5) The temperature increment due to adding heat 
from fuel combustion between 6 and 10 is determined. 

(6) Gas state conditions at 10 are calculated. 

(7) The gross thrust from the overall configuration 
is found. 


From a simplified physical point of view, for nozzle 
exit conditions below or slightly above sonic, the in- 
crease in gross thrust developed by the afterburner over 
that of the engine alone may be interpreted as being due 
to an increase in nozzle exit area over which the jet 
total pressure acts, since the difference between Pio’ 
and Py is slight and Ay is much larger than Ay’. 
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Method of Analysis 


In this section the complete method of analysis is 
described in detail, including all of the necessary for- 
mulas. Readers not acquainted with some forms pre- 
sented are referred to the appendix for development of 
the less familiar equations. 

Since the conditions at station 5 for a given ram ratio, 
air speed, and altitude are to be constant and equal to 
those at 5’ for the same engine operating point, the 
first step is the calculation of these conditions from the 
known engine performance data, so that they may be 
used in the actual afterburner analysis. That is, the 
total pressure and temperature at the nozzle exit of the 
standard tailpipe are computed, and assumed or calcu- 
lated tailpipe total pressure losses are subtracted from 
the former to obtain P;’. The total pressure loss be- 
tween stations 5’ and 10’ in tailpipes less than 5 ft. 
long varies between 3 per cent and 5 per cent. T;’ is 
assumed equal to Tyo’ generally, since little total heat 
content is lost between the two stations; a more rigor- 
ous solution may include the temperature decrement 
from radiation and convection. 

Before Pio’ can be found for the exit nozzle, it is nec- 
essary to determine whether the flow at the exit is sub- 
sonic or sonic, since the latter implies a supercritical 
value of the pressure ratio Pjo'’/p,, requiring different 
pressure formulas. 


Mj = VFj'/7pAw’ (1) 


where M,’ is the theoretical Mach Number attainable 
by the jet, corresponding to the pressure ratio Pyo’/p,. 
This equation is valid only for Mach Numbers less 
than or equal to 1.0; values of M,’ greater than unity 
are not true values but do indicate sonic nozzle throat 
velocity, which is actually the only information de- 
sired. 
If < 18: 


Pro! = pol [Fy'(y — + (2) 
If M,’ > 1.0: 
+ [(y - 
= (45 +5 (3 
Then P;’ may be found: 
Ps! = 100P//(109 — % loss) (4) 


Now, at any point in the gas stream, the Mach Num- 
ber, gas flow area, total pressure, and total temperature 
are related: 


W,VT MV 6) 
where 


W, = Wa + 


A curve of a vs. M is given in Fig. 2. Variation of 
this curve, and the subsequent one for 8, with tempera- 
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MACH NUMBER, M 


Fic. 2. Gas flow parameters. 


ture due to change in y is very slight; separate curves 
may be drawn for several temperature ranges, but the 
error in neglecting this is small. As a general rule, the 
value of ai’ is not compatible with the function of 
My’ in Eq. (5) using engine performance data. That 
is, substitution of the performance values of gas flow, 
total temperature, nozzle exit area, and total pressure 
at station 10’ do not give a value of ay’ (and thus Mjp’ 
from Fig. 2) equal to that obtained by calculating 
My’ from Eq. (1). W,’, P10’, Aw’, and F;’ are assumed 
to be reliable, and, since M,’ is a function of them, the 
values of Ti’ (or 7;’) given in the performance data for 
the engine are not used. Tyo’ is calculated instead by 
finding a’ from the graph at Mio’, where Mi’ = M,’ 
up to M,’ = 1.0, and retains a value of unity for higher 
values of M,’ [calculated from Eq. (1)]. 


Then, 
P9'A10'\? 

At this point the standard tailpipe and fixed nozzle 
may be replaced with the afterburner and variable 
nozzle. The primes are removed for the afterburner 
quantities. 

As indicated, the basic premise for the afterburner 
operation is that all conditiens at station 5 remain 
equal to those at 5’ at a given engine speed, airplane 
speed, and altitude. Or 


P; = Py’, Ts = Ty’, 


Friction and diffusion losses in total head must be 


assumed or calculated between stations 5 and 6. Then, 
P., = P;(100 — % loss)/100 (7) 

and 
as = W,VTs/PoAs (8) 


Mg is then found from Fig. 2. Now, 
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M 
VygRT V1 [(y — 1/2]? 


This curve is also presented in Fig. 2. From the 
value of Ms, Bs may be found. Then, 


Ve = BeV ygRT (10) 


The foregoing determination of velocity at station 6 
is not necessary for the actual thrust solution, but it is 
essential in practical design problems to ensure having 
V¢ low enough for maintenance of combustion—that is, 
to avoid flame blowout. 

The major portion of the engine operating time will 
be with the afterburner unlit, yet causing a total pres. 
sure loss in the jet. The thrust available at this con- 
dition may be found as follows: 

The total pressure loss between stations 6 and 10 
must be determined. Then, 


Py = P.(100 — % loss)/100 (11) 


As before, 1, must be found, though by a different 
formula, since Pi) is known instead of F;. Thus, 


M, = V [2/(y — — 1] (12) 


Reference is made to the a vs. M curve, following the 
practice previously outlined for M, > 1.0, and ay 
is found. Then, 


An = WV To/ P. (13) 
The jet thrust becomes: For M; < 1.0, 
F,= [(Pro/ — 1)[2vAwp./(y — 1)] (14) 
for 1.0, 


(9) 


+ 1) 
+ [(y — 


If the afterburner is lighted the following analysis 
is appended to the conclusion after Eq. 10: 

The fuel flow may be selected in any arbitrary man- 
ner up to a maximum value of 


= (0.067W. — Wrz.3.)/n (16) 


This maximum value assumes that all of the excess 
air from the turbojet exhaust is burned. Then the 
temperature rise from station 6 to 10 is 


AT = + Wrro.) (17) 


Eq. (17) is, of course, an approximation to the actual 
integrated temperature rise during the combustion 
process. 

Then, 


p.) Ay (15) 


Ti = Ts + AT (18) 


The momentum of combustion loss and that due to 
friction must be considered between stations 6 and 10; 
then, 
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Pi = P.(100 — % loss) /100 (19) 


As before, M, is found by Eq. (12), using y values 
specified approximately by Ty (accurately by tp). 
Then a1 is taken from Fig. 2, using the aforementioned 
rule for M, > 1.0, and 


Aw = Tw/Pwow 


The jet thrust is found from Eqs. (14) and (15), 
using appropriate values of y. Net thrusts in all calcu- 
lations are found, when desired, by 


Of course, the general method outlined may be al- 
tered to suit a particular problem. Often the combus- 
tion-chamber area is to be determined at a given condi- 
tion for a maximum value of Vs, compatible with steady 
combustion. In such a case, Eqs. (8), (9), and (10) 
are reversed in order, and solution is made for Ag. 
Once selected, As must remain constant for practical 
analysis of a given unit. The afterburner fuel flow 
may be established by any of several criteria. Possible 
ones include: (1) constant AT for all conditions, (2) 
constant Ty for all conditions and (3) constant after- 
burner fuel-air ratio. 


(20) 


(21) 


EXAMPLE APPLICATION 


Since the formulas used in this investigation are 
straightforward thermodynamic derivations, the author 
feels that extensive sample calculations would be super- 
fluous. However, in order to emphasize the enormous 
benefits attainable in certain phases of aircraft perfor- 
mance with an afterburner installation, an analysis has 
been made of a hypothetical airplane having the follow- 
ing characteristics: 


Wing area 300 sq.ft. 
Aspect ratio 5.09 
Span 39.1 ft. 
Sweep of quarter-chord line 30° 
Minimum parasite drag coeffi- 

cient 0.018 
Normal gross weight 18,000 Ibs. 


The airplane is powered by four hypothetical tur- 
bojets, supplemented by afterburners operating with 
the following assumed characteristics: » = 75 per 
cent, Tio = 4,000°R, and Ayo varied to maintain condi- 
tions at station 5 equal to those at 5’ in the basic 
analysis. 

Fig. 3 shows the thrust required curve at sea level as 
a percentage of standard sea-level static thrust at 
military (take-off) r.p.m. plotted against true air speed 
as a percentage of maximum speed at sea level at mili- 
tary power. The thrust available curves, with and 


without afterburning, are plotted in similar percent- 
ages. It is apparent that the maximum speed at sea 


level without the afterburner occurs at a Mach Number 
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Fic. 3. Sea-level thrusts required and available. 


slightly above that of drag divergence for the airplane. 
Thus the increases in maximum speed when the after- 
burner is lighted is only 10 per cent above the turbojet- 
alone top speed; this is still a considerable increment 
in the transonic zone. The variation of maximum 
speed with altitude appears, in percentages, in Fig. 4. 

The manifold increase in excess thrust available for 
climb is apparent in Fig. 3; this increment, combined 
with the increased speed for best climb, is reflected in 
the rate-of-climb versus altitude curve of Fig. 4, which 
indicates a 200 per cent increase in sea-level rate of 
climb. Similar gains in one- or two-engined-inopera- 
tive climb performance accrue, introducing a note- 
worthy performance safety factor for pure jet aircraft. 

But probably the most important benefit derived 
from the afterburner is indicated by a study of the dis- 
tances to take off and climb 50 ft., the results of which 
are shown in Table 1. The first two distances, marked 
“normal take-off distance,’’ were calculated by conven- 
tional means, including the Welling transition equa- 


tions. The second two distances, labeled “emergency 
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Fic. 4. Performance summary chart. 
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take-off distance,” were calculated by the same miethod, 
but one engine was assumed to fail at the point where 
the airplane reached take-off velocity, and the transi- 
tion and climb to 50 ft. was conducted with only three 
engines (or engines and afterburners) functioning. 
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TABLE 1 
Turbojets at 
Turbojets Alone Military 
at Military Power with 
Power Afterburners 
Normal take-off 
distance, ft. 4,520 2,255 
Emergency take-off 
distance, ft. 8,650 2,825 


A wing loading of 60 Ibs. per sq.ft. was used, corre- 
sponding to a fully loaded take-off condition. It is ap- 
parent that the difference between the afterburner on- 
and-off configurations represents the difference be- 
tween a safe or a dubious take-off, particularly if a 
commercial transport is being considered. In the ex- 
ample case, three engines with afterburners operating 
have about 15 per cent more thrust than four engines 
without augmentation. 


Appendix 


DERIVATION OF FORMULAS 


When the pressure ratio P1o/p, is subcritical, a vena 
contracta exists somewhat downstream of the nozzle 
exit; at the contracta, and there only, the static pres- 
sure in the jet is equal to ambient static pressure, p,. 
Where this is true, the area-pressure term may be omit- 
ted from the jet-thrust expression (since the pressures 
are balanced), and the jet thrust becomes: 


= 


This may be transformed into a more usable form as 
follows: since 


(A-1) 


ms = 
and 
ay = = YPo/p; 
Then, 
Fy = pV PA; 

= p;A;M/a/ (A-2) 

= 
Also, since 


Pi/bo = {1 + [(y — 
Then, 
M? = [2/(y — — 


In practice the variation in jet static pressure and 
area between the nozzle exit and the vena contracta 


1} (A-3) 


has been found to be slight, on the order of 2 to 5 per 
cent. And since in this analysis the assumption of no 
vena contracta is applied to both the turbojet alone 
and the engine with an afterburner, the error of assump. 
tion is compensatory. 

Assuming, then, that A; and M; are equal to A, and 
My, respectively, Eq. (A-3) may be substituted in Eg, 
(A-2) and rearranged: 


[Fi(y — 1)/2ypAwl + 1 = 
or 
Pro = pol [Fi(y — 1)/2vAwpo] + 
If Pio is known and F, is wanted, then 
= — — 1)] (A) 


It should be noted that, though Eqs. (A-2) and 
(A-3) were developed on the basis of a subcritical pres. 
sure ratio, the latter is accurate for supercritical ratios 
as well. Since it is also necessary in some instances to 
determine whether or not Vj is sonic, knowing only F, 
and A,, Eq. (A-2) may be solved for M;, and the value 
obtained is accurate up to M,; = 1.0. Values of M, 
greater than unity are not true values but do indicate 
that the nozzle velocity is sonic. This is actually the 
only information desired, since it indicates whether 
Eq. (A-4) or Eq. (A-5) may be used (when /,; < 1.0) or 
if other equations, developed subsequently, for super- 
critical conditions must be applied. 

Now, if M,; > 1.0, the following basic equation must 
be applied at station 10, which includes the pressure- 
area term to account for the difference in static pres- 
sure of the jet and ambient conditions: 


F; = (W,/g) Vi + — Po) (A-6) 
But Vio = dy, since My = 1.0. Therefore, 
W, 
F; = + — fo) (A-7) 
Now 
W, = pwA 10 Viog 
= (ygbhwA 
But My = 1.0. Therefore, 
Pio = Waw/ygAw (A-8) 


Substituting this in Eq. (A-7), 
F; = (Weaw/g){1 + (1/7)] — Apo 
Solving for a, 
[(F; + + 1)] 
which may be substituted in Eq. (A-8) 
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Pyo/pro = {1 + [Cy — 1)/2] P 
Therefore, 


+ PoAw) [1 + (¥ 


Pw 


Ay y¥+1 (A-9) 
( Fs p.) 
Ax y¥+1 
Or, if Pio is known and F; is desired 
(A-10) 


In reference 1, Bailey presents the formula 


TURBOJET ENGINES 311 
then, 
WT _ 
- 
(A-12) 
From the same reference, 
M = V/V ygRt (A-13) 
Since 
t = T/{1 + [(y — 1)/2]M} 
then, 
V M 


VygRT) V1 +4 [(y — 1/2)? 


In the September, 1947, issue of the JouRNAL, William F. Milli- 
ken, Jr., presented an excellent paper on ‘‘Progress in Dynamic 
Stability and Control Research.”” The importance of carefully 
planned and instrumented flight-testing programs to the success- 
ful development of airplanes and missiles was emphasized. 

Milliken showed that the changes in forward speed were small 
enough, for the frequency range he was considering, to enable him 
toneglect the X-axis equation of motion in the analysis. Further- 
more, in order to simplify the analysis, he assumed the derivatives 
Z,, Mix, and Zs were zero. Although each of these derivatives is 
small, unfortunately there are cases in which they combine in such 
amanner that their total effect is far from negligible. For ex- 


= M (1 + we) (A-11) REFERENCE 
Since 1 Bailey, Neil P., The Thermodynamics of Air at High Velocities, 
/a=1) Journal of the Aeronautical Sciences, Vol. 2, No. 3, p. 227, July, 
pb = P/{1 + [(y — 1944. 
Letter to the Editor 
Dear Sir: ample, suppose we include these three additional derivatives in 


the equations of motion used by Milliken. Using the same 
method of analysis he used, we obtain an expression for the ratio 
of the maximum amplitudes of the normal acceleration || to the 
maximum amplitudes of the pitching velocity |g|. This ratio 
|n|/|q| was used by Milliken to determine the slope of the lift 
curve Z, by the expression (on page 506): 


|n| UZ. U 
lal +02 V14 (w/Ze)* 


(1) 


Using the same dimensional notation as in the subject report, 
we obtain for |n|/|q|, 


In| _ bi wb (ki — (2) 
k 2\2 
lq! — Cus’)? + ] 
wh; 
where Note that if Z3, Z,, and Mj; all equal zero, 
b = —(M,+ + A=a=G= (3) 
k = UMy Aa=G=-C=0 
by = (b/G) — (MuaZq/Cr) and Eq. (2) reduces to Eq. (1). 
ky = (k/G) — (MwZ,/Ci) Obviously Eq. (1) is much simpler and easier to use than Eq. 
CG =1— Mi(U + Z,) (2). It is our present purpose to estimate values of the con- 
C, = 1 — (2,Mu/MsZw) stants C; to C, and and in order to evaluate the accumulated 
C; = 1 + (25Mie/MoZw) errors of making the assumptions of Eq. (3). Unfortunately, 
Cy = Z5Mi/MsZw Zw, Mw, as well as other derivatives we should like to be able to 
ZsUMis _ obtain explicitly, are contained in these constants. If, however, 
ay a nag Ms Ms we consider the constants C; to C; as small correction terms, then 
= 1 25 1 — UM3) they can be estimated satisfactorily either theoretically or from 
ae: UZ ( wind-tunnel data. 


| 
: 
| 


312 JOURNAL OF THE AERONAUTICAL SCIENCES—MAY, 1948 


T 
T 
OF E SJECT RE T(Fig 14) 
€.6. (MAC), OF 22%, 28%8% 
Lit f 
| } 
to 9 z 5 | 
TH 


FIGURE |: EFFECT OF 25,24, &M;, ON THE VARIATION OF 2,, WITH 
FREQUENCY OF OSCILLATION. B-25u AIRPLANE 175M.P.H. EQUIVALENT 
AIR SPEED. 10,000-FT. PRESSURE ALTITUDE. 


Let us evaluate Eq. (2) for the conditions of Fig. 14 of the sub- 
ject report: equivalent air speed of 175 m.p.h., 10,000-ft. pres- 
sure altitude, and various c.g. locations from 22 to 30 per cent 
of the M.A.C. Most of the values of the derivatives are obtained 
directly from the subject report. The derivatives M;, and Mz 
are derived from the rather accurate theory of N.A.C.A. T.N. 
No. 771 (Robert T. Jones and Leo F. Fehlner, ‘“‘Transient Effects 
of the Wing Wake on the Horizontal Tail’). It might be noted 
that the results of 7.N. No. 771 yield much larger values of ef- 
fective time lag than is usually considered. (At usually equals 
tail length divided by velocity.) 

The values of Z, obtained using Eq. (2) are compared in Fig. 1 
with those of Eq. (1) (Fig. 14 of subject report). It is clear that 
the values of Z,, obtained using Eq. (2) agree more closely with 
the wind-tunnel values than the values of Z from Eq. (1). How- 


ever, the scatter of points is unchanged. This indicates a lack of 
accuracy of the instrumentation, since Z» certainly should not be 
a function of c.g. position. 

The effects of Z, and Mj on the estimated values of Z~ proved 
to be negligible for this airplane configuration and frequency 
range. For other airplane configurations or speed ranges, or for 
evaluating other derivatives, Z, and Mj correction terms age 
important. For example, 7; changes the calculated values of 
damping and spring constants, }; and k,, by about 6 per cent, 
Analysis of various modern dynamic problems, such as the aboyg 
has shown even fairly small terms such as Z;, when multiplied by 
w or w? (if = |6| sin wt, = —w?|5| sin wt) may be important fop 
large values of w. The importance of the Z; term is emphasize 
in that the apparent values of Z,, from Eq. (1) increased with fre 
quency, whereas the inclusion of the Z; term resulted in Z, aps 
parently decreasing with frequency. The physical reason for the 
importance of the Zs term at the higher values of the oscillatiog 
frequency is that, at these high frequencies, the aircraft does not 
have time to change its angle of attack very much during the shor 
period of such acycle. Therefore, practically the total normal ag 
celeration |”| results from the elevator lift rather than from the 
wing lift. An erroneous answer for |m|, and consequently Zy 
therefore results if Zs; is not considered at the higher oscillation 
frequencies. (Recent conversation with Mr. Milliken indicates 
that after publication of the subject report he had independently 
discovered the importance of the Z; term.) The only safé rule 
appears to be to include all possible terms and discard them only 
in the final equations if they prove small enough. This procedurg 
unfortunately involves an excessive amount of work in most caseg, 
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